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Summary 


The energies of some neutron groups leading to low excited states of “C in the 
reaction B(d,ny)"C have been measured. These lead to values of 4-3--0-3 MeV and 
6-53+-0-02 MeV for the energies of the second and fourth excited states respectively. 
y-Ray spectra have also been studied in coincidence with different neutron groups. 
Ground state transitions were observed from the second, third, and fourth excited 
states, together with cascade decays of the fourth excited state through each of the 
second and third. Deductions from these y-ray spectra considerably reduce the number 
of spin possibilities which have been found for these levels by other workers. 


I. INTRODUCTION 
A good deal of information about the energy levels of a nucleus can be 
obtained from the observation of y-ray transitions between such levels. In 
the case of UC, the parities of the first few levels are known, and measurements 
of the relative intensities of the possible y-ray transitions between these levels 
make it possible to deduce the most probable multipolarities of the y-rays and 
hence to select the most probable spin values from the many alternatives which 

have been provided from other measurements. 


The current knowledge regarding the energies, spins, and parities of the level§ 
of 11C is listed fully in the compilation of Ajzenberg-Selove and Lauritsen (1959), 
and the information relevant to the work described in this paper is reproduced, 
together with the equivalent data on the mirror nucleus “B, in Figure 1. The 
parities of all these levels are known to be odd from the /,=1 stripping patterns 
found by Cerineo (1956) and Maslin, Calvert, and Jaffe (1956) for the neutron 
groups in the reaction B(d,n)"C. These stripping data also restrict the spins 
of all these levels to 3/2, 5/2, 7/2, or 9/2, and in the case of the ground state, the 
6+ decay evidence restricts the spin to 1/2, 3/2, or 5/2. The value 3/2 is con- 
sidered the most likely, this being the spin of the “B ground state. 

Apart from the 2 MeV y-ray from the first excited state, the only y-rays 
which have been observed from C are those observed by Bent et al. (1955) and 
Sample ef al. (1955) in the presence of all other y-rays following the deuteron 
bombardment of B, and the only certain assignment is the ground state 
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transition from the 6-50 MeV state. In the work to be described below, the 
spectra of y-rays involved in the decay of individual levels are examined 
separately, their observation being in coincidence with the emission of selected 
neutron groups in the B(d,ny)"C reaction. 

The situation in the mirror nucleus, 1B, where the ground state spin is known 
to be 3/2-, is a little clearer. The /,=1 stripping patterns found by Evans and 
Parkinson (1954) in the reaction B(d,p)"B place the same restrictions on the 
spins and parities of the low levels of “B as the /,=1 patterns did in the case of 
110, but these spin possibilities have been considerably reduced in most cases by 
the angular distribution and angular correlation measurements of Thirion (1953), 


6-50 3/27 To of2- 
“75 


Fig. 1.—Accumulated data concerning the energies and spins 
of excited states of “C and ™B. 


Bair, Kington, and Willard (1955), and Ferguson e¢ al. (1958). Of particular 
interest is the conclusion that the spin of the 2-13 MeV level is 1/2, in agreement 
with the intermediate coupling model, but in contradiction to the 1,=1 stripping 
data, the contradiction being resolved by Wilkinson (1957) by postulating spin 
flip stripping. Since the situation in the case of the 2-01 MeV level of 1C is exactly 
analogous, it would be reasonable to consider spin 1/2 as a possibility for this 
level also. 

Our experiment concerns only “OC, consisting of measurements on the 
reaction “B(d,ny)"C. Two types of measurement have been made: firstly 
the energies of the neutron groups were measured by a time-of-flight motuons 
and secondly, the y-ray spectra in coincidence with the individual neutron eae 
were observed. 
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II. EXPERIMENTAL DETAILS 

A beam of deuterons from the 750 kV electrostatic generator struck a 
separated boron target (99 per cent. °B, supplied by A.E.R.E., Harwell) deposited 
on a nickel backing. Two target thicknesses were used: 90 ug cm for the 
neutron energy measurements and 260 ug cm for the y-ray measurements. 
y-Rays were detected by a NalI(T1) crystal placed as close as possible to the target 
and at 90° to the direction of the deuteron beam. This crystal was coupled to 
an RCA 6342 photomultiplier tube. Two crystal sizes, 14 in. diameter by 1 in. 
and 3in. diameter by 3in. were used. The neutron detector consisted of a 
3in. diameter by ? in. thick cylinder of “ Pamelon ”’ plastic scintillator (manu- 
factured by Isotope Developments Ltd.) coupled via a “ Perspex ” light pipe 
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Fig. 2.—Block diagram of delayed coincidence circuit. 
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to an EMI 6262 photomultiplier. This counter was surrounded by a lead and 
paraffin wax shield with a collimator defining the directions of the fast neutrons 
incident on the scintillator. The whole of this assembly could be moved so 
that both the length and direction of the neutron flight path were easily varied. 


While the neutron counter was always located in the horizontal plane 
containing the incident beam direction, different locations were used for the two 
Nal(T1) crystals. The 14 by 1 in. crystal was placed with its axis in this same 
horizontal plane but, when the 3 by 3 in. crystal was used, it was located below 
the target with its axis vertical. 

The block diagram (Fig. 2) shows the time-of-flight circuitry. By using the 
largest recommended value of the high voltage applied to the 6342 photomultiplier 
and a small (1200 Q) collector load, a pulse from the y-counter could be obtained 
which was fast enough to operate the fast coincidence unit (Harwell type 1153A) 
at its shortest resolving time of 5mys. With this arrangement, the delay 


‘ 
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distribution used for testing—that obtained with the counters close together and 
a source of 22Na between them—had a width at half height of 7 mus. 

The interpretation of a delay distribution for a given neutron flight path 
requires a knowledge of the location, on the variable delay axis, of the ‘ zero 
delay ” position, i.e. the setting of the variable delay for which pulses from events 
in the two scintillators that are truly prompt are brought into coincidence. In 
our case, because of the severe differentiation of the collector pulses from the 
y-counter, the rise time of these pulses varied with their magnitude. This 
resulted in a small variation in the inherent delay of such pulses, so that the 
zero delay position was dependent on the energies of the y-rays involved in the 
coincidences. For this reason, it was not possible to use the position of the peak 
observed in the test delay distributions taken with a ?*Na source as the zero delay 
position for the (n,y) coincidences, nor was it possible to use the position of the 
prompt peak in an actual run, this latter peak being contributed to by coincidences 
such as y-y cascades following B(d,p)"4B* which involve a large range of y 
energies, and ‘“‘ inverted’ (n,y) coincidences in which the neutron is detected 
in the y-counter and vice versa. Hence the measurements of neutron energy 
were obtained by finding the position of the appropriate delay peak at various 
lengths of flight path and determining the neutron velocity from the slope of the 
resulting distance-time graph, all points on which were then produced by 
coincidences in which the same y-ray was involved and were measured with 
identical circuitry. The discriminator settings in the slow coincidence circuit 
were chosen to eliminate, as far as possible, events due to coincidences, both 
random and true, other than those under observation. 


IIT. DELAY DISTRIBUTIONS 
(a) The Neutron Group to the 6-50 MeV State of UC 

Delay distributions were measured at neutron directions of 0 and 90° relative 
to that of the incident deuteron beam and at effective deuteron energies, allowing 
for a mean energy loss of 20 keV in the target, of 580 and 630 keV. The neutrons 
of this group are of low energy (about 0-5 MeV) so that the pulses produced by 
them in the neutron detector were of the same order of size as the noise pulses. 
Thus the lowest available setting of the discriminator in the neutron channel 
had to be used. In the y-channel, the discriminator was set to exclude pulses 
produced by electrons of energy up to 2 MeV from the 1°8I activity in the Nal 
(see Section IV (a)), thereby favouring the detection of coincidences associated 
with the relatively intense 6-5 MeV radiation (Bent et al. 1955). In this way a 
true-to-random ratio of about 1 was obtained. A typical delay distribution is 
shown in Figure 3. Figure 4 shows the delay versus distance relationship for 
one set of delay distributions, and the results obtained from these measurements 
are shown in Table 1. 


(b) The Neutron Group to the 4-26 MeV State of 40 
For the 0° direction and a deuteron energy of 580 keV the neutron groups 
to the 2-01, 4-26, and 4-75 MeV states of C should have energies of approxi- 
mately 4:9, 2-8, and 2-3 MeV respectively, with the corresponding flight times 
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Les metre of 33, 44, and 48 mus. With delay distributions having a width at 
alf height of 10 mus, Tauber long flight paths would be required to resolve these 
groups. Under these conditions the count rates in our detector would be pro- 
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Fig. 3.—-Typical delay distributions. Curve (i) neutron group to 
6-50 MeV state (90° direction; flight path 33-:6cm). Curve (ii) 
neutron group to 4:26 MeV state (0° direction ; flight path 50-8 cm). 
The random rate was measured by inserting a large delay in the neutron 
channel, and the higher count rate at large y-channel delays is 
attributed to true coincidences involving scattered neutrons. 
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Fig. 4.—Delay v. flight path measurement for neutron group to the 6-50 MeV 
state (neutron direction 90°; beam energy 650 keV). 


hibitively small for a series of delay distributions of the type hitherto described. 
However, the neutron yield measurements of Paris and Endt (1954) indicate 


that at this energy the intensities of the groups in the 0° direction should be in 
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TABLE 1 
RESULTS OF DELAY DISTRIBUTION MEASUREMENTS FOR THE SECOND AND FOURTH EXCITED 
STATES OF !'C 


Slope of 
Deuteron Neutron Delay- Q Energy of 
Energy Counter Distance (MeV) Level 
(keV) Angle Curve (MeV) 
(mus m-) 
580 0° 4442 2-16+0-27 4-3 +0-3 
580 0° 104+4 —0-08-+0-04 6-55+0-04 
580 90° 115+5 —0-04+0-04 6-51+0-04 
630 0° 98+2 —0:07+0-02 6-54+0-02 


630 90° 111+3 —0-06+0-03 6-53+0-03 


Mean Value: 6-53+0-02 


the ratio 6:1:0-3, the 2-8 MeV group being the strongest and the 4:9 MeV 
group the weakest. Thus, a set of delay distributions taken under these con- 
ditions can be interpreted as giving a reasonable value of the energy of the 
4-26 MeV state only. 
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Fig. 5.—Delay distribution of neutrons to the 4-26 and 4-75 MeV 
states (neutron direction 90°; flight path 171-5 cm). 


In this case, the discriminator setting in the neutron channel could 
be increased to remove noise pulses and also those due to the annihilation radiation 
of the “C. The resulting improvement in the true-to-random ratio and in 
resolution is evident in Figure 3. Table 1 includes the results of the Single set 
of delay distributions made on this group. 
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(c) The Neutron Group to the 4-75 MeV State of UC 

In the interpretation of the coincidence y-ray spectra (Section IV (c)) it is 
desirable to have an estimate of the relative intensities of the neutron groups to 
the 4-26 and 4-75 MeV states. To provide such an estimate, a single delay 
distribution was taken over a flight path of 171-5 cm in the 90° direction. The 
angular distribution measurements of Graue (1956) indicate that in this direction 
the intensity of the weaker group (to the 4-75 MeV State) should be a maximum, 
while that of the stronger group should be at its minimum value. The result 
of this measurement (a 30-hr run), performed with a 260 ug cm target and a 
beam current of 0-5 uA at 600 keV bombarding energy, is shown in Figure 5. 
Comparison with the delay distribution shown in Figure 3 for the (effectively) 
monoenergetic neutron group of approximately the same energy shows the 
presence of two unresolved neutron groups. The separation of their two delay 
peaks is estimated to be 9+2 mus, corresponding to an energy difference of 
0-5+0-1 MeV, in agreement with other determinations of the energies of these 
two levels, and the ratio of their intensities is approximately 2, the group going 
to the second excited state being the more intense. 


IV. COINCIDENCE y-RAY SPECTRA 
(a) Energy Calibration 

The spectra of the y-rays in coincidence with the various neutron groups 
were recorded in a 100 channel pulse-height analyser (Sunvic type PHA2). 
Calibration points were obtained from the spectra of the 2-62 MeV y-ray of 
208Pb (from a radiothorium source) and the 4-43 MeV y-ray of ##C (from a Po-Be 
source). With the beam striking the target the count rate in the y-counter 
was very high, with large contributions from the “C (20-min half-life) annihilation 
radiation from the target and the neutron-induced ?°I (25-min half-life) activity 
in the Nal(T1) crystal itself. This high count rate, estimated to be in excess of 
105 per sec, had to be tolerated in order to obtain usable coincidence count rates. 
The resultant pile-up of pulses was sufficient to make the pulse-height resolution 
of the counter a good deal worse than could be obtained at low count rates, and 
also slightly reduced the gain of the photomultiplier. For this reason, calibration 
spectra were taken several times during each run, the beam being removed from 
the target only for the few minutes required to insert the radioactive sources and 
record the spectra, so that, under these conditions, the total count rate in the 
counter was still very high. Typical spectra are shown in Figure 6. 

In addition, the energy of the ground state transition from the 6-50 MeV 
state could be regarded as being known from our delay distribution measurements 
to a better accuracy than we could hope to achieve in our scintillation spectra. 
Hence, in the coincidence spectra in which this y-ray appeared, gon a 
further calibration point, obtained during the actual running conditions for the 
coincidence spectrum. ‘The fact that this point, together with the two obtained 
from radiothorium and Po-Be, gave a pulse-height versus energy relation which 
was linear within the experimental accuracy was rogarded as a check on the 
validity of the calibration procedure using the radioactive sources described 


above. 
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(b) y-Rays from the 4-26 MeV Level of UC 
With the neutron counter in the 0° direction and a short flight path, the 
y-rays in coincidence with the unresolved neutron groups to the 4-26 and 
4-75 MeV states will be principally those due to de-excitation of the second 
excited state. A spectrum obtained using the small NalI(T1) crystal is shown in 
Figure 7. Comparison of this figure with Figure 6 (b) clearly suggests the presence 
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Fig. 6.—Calibration y-ray spectra. (a) RdTh, 1} by 1 in. crystal, (b) Po-Be, 14 by 1 in. 
crystal, (c) RdTh, 3 by 3in. crystal, (d) Po-Be, 3 by 3in. crystal. 


of only one y-ray. The poorer resolution probably results from a combination 
of the higher count rate and small random drifts in gain during the run. For the 
runs reported here, no systematic drifts in the calibration spectra were observed, 
although random drifts of up to +-1 channel did occur at times. 


The main peak in Figure 7 is interpreted as the two-quantum escape peak 
of a y-ray of energy 4-2+0-1 MeV, and this is identified as the ground state 
transition from the 4:26 MeV state. A cascade from this level through the 
2-01 MeV state would produce two y-rays of energy in the vicinity of 2 MeV. 
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From Figure 7 we conclude that the intensity of any 2 MeV radiation which 
may be present is less than 10 per cent. of that of the 4-2 MeV y-ray and hence 
the transition probability for de-excitation via the 2-01 MeV state ig less than 
5 per cent. of that for the ground state transition. However, it will be seen, in 
the following section, that this upper limit may be reduced to 24 per cent. 


(c) y-Rays from the 4-75 MeV State Ope LG 
When the neutron counter is in the 90° position, our delay distribution 
measurements have shown that approximately one-third of the true coincidences 
in the composite delay peak of the neutron groups to the 4:26 and 4-75 MeV 
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Fig. 7.—y-Ray spectrum from 4-26 MeV state (small NaI(T1) 
crystal), 


states are associated with the 4-75 MeV state. Thus the y-ray spectrum in 
coincidence with this composite delay peak may provide some information 
about the de-excitation of this level. 

With a short flight path and using the small Nal(T1) crystal the spectrum 
shown in Figure 8 was obtained. It is not possible to select either of the two 
curves shown in this figure as being a better fit to the experimental points than 
the other and it is clear that very good statistics indeed would be needed to be 
able to decide between them. Thus, from this measurement we are not able to 
discover whether or not the 4-75 MeV state is de-excited to any extent by a 

ion. 
. Stam ne ana under the same conditions but using the large NalI(T1) 
crystal is shown in Figure 9. The “ smearing out ” of the spectrum from a 
crystal of this size by the use of a high count rate is clearly seen by reference to 
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Figure 6 (d). It is also apparent from Figure 7 that this effect is less serious in 
the spectrum from the small crystal, partly due to the lower count rate and 
partly due to the nature of the latter spectrum in which the two-quantum escape 
peak is dominant at these energies. 
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Fig. 8.—y-Ray spectrum from 4:26 and 4:75 MeV states 

(small NaI(Tl) crystal). The full curve is the spectrum expected 

if one-third of the y-rays were 4-75 MeV and two-thirds 

4-26 MeV. The dotted curve is the spectrum for 4-26 MeV 
alone. 


Figure 9 again provides no evidence for or against the presence of a 4-8 MeV 
y-ray but does confirm the low intensity of any cascade through the 2-01 MeV 
state. One could perhaps suggest the presence of a peak having a height of the 
order of 5 per cent. of that of the main peak in the vicinity of 2 MeV, as shown 
by the dotted curve in Figure 9, but it would hardly be justified by the statistics 
of the points in this region. Thus it seems reasonable to set the upper limit for 
de-excitation of the 4-26 and 4:75 MeV states combined via the 2-01 MeV state 
as 23 per cent. of the total probability of de-excitation of these levels either by 
direct ground state transitions or by cascade from the 4-75 MeV to the 4-26 MeV 
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state and thence to the ground state. Hence less than about 74 per cent. of the 
transitions from the 4-75 MeV excited state alone are from this level direct to 
the first excited state. However, the measurements provide no information 
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Fig. 9.—y-Ray spectrum from 4-26 and 4-75 MeV states (large 
Nal(T1) crystal). 


about the probability of a cascade from the 4:75 MeV to the 4:26 MeV state. 
It was not possible to obtain any observations of neutrons in coincidence with a 
0-5-MeV y-ray (the first transition in such a cascade) because of the high 
annihilation radiation background. 

TABLE 2 


PEAKS IN THE SPECTRA OF Y-RAYS IN COINCIDENCE WITH THE NEUTRON GROUP TO 
THE FOURTH EXCITED STATE OF UC (smauL Nal(Tl) cRySTAL MEASUREMENTS) 


Peak Energy 
(MeV) Identification 
Run 1 Run 2 Mean 

* ney G9) 

Y 2-3+0-1 2-3 2-3 MeV; full energy peak 
3°2+0-1 3-4+0:1 3°3 4-3 MeV; double escape peak 

f4:3 MeV; single escape peak 
Sadie Bare Oct A ‘4:8 MeV; double escape peak 
4-340:2 4-3-+0-15 4-3 4:8 MeV; single escape peakt 
6-50 6-50 6:50 MeV; double escape peak 
(calibration point) 


* Spectrum unreliable in this region due to electronic fault. 
+ Plus a small contribution from the 4:3 MeV full energy peak. 


(d) y-Rays from the 6 -50-MeV State 
The spectrum of y-rays in coincidence with the low energy neutron group 
was observed with both the small and large Nal(T1) crystals. The small crystal 
results are shown in Table 2, and one of the small crystal spectra is shown in 
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Figure 10. It will be seen that there is clear evidence for a cascade of two ytaye 
of approximate energies 2-3 and 4-3 MeV. Using the best available level energies 
as shown in Figure 1, the two possible cascades which might fit this description 


ay (i) 4-49 MeV followed by 2-01 MeV, through the first excited state, and 


(ii) 2-24 MeV followed by 4-26 MeV, through the second excited state. 


We believe that the accuracy of our y-ray measurements is sufficient to 
exclude the first of these possibilities. The peak falling at 2-3 MeV in Figure 10 
is well defined and appears at the same energy in the coincidence spectrum 
obtained with the large NaI(T1) crystal (Fig. 11). Our energy calibration seems 
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Fig. 10.—y-Ray spectrum from 6:50 MeV state (small NaI(T1) crystal). 


reliable enough to exclude the possibility that this peak could be due to a 2-0 MeV 
y-ray. In the same way, the peaks in the middle energy range cannot be fitted 
to a 4:5 MeV y-ray, particularly in view of the proximity of our 4-43 MeV 
calibration points. Thus it seems that the predominating cascade from the fourth 
excited state is via the second excited state, and that the probability of occurrence 
of this cascade is somewhat less than, but of the same order of magnitude as 
that of the direct ground state transition. 


There is also evidence, more particularly from the small crystal spectra, 
for the presence of a 4-8 MeV y-ray, corresponding to the 4-75 MeV ground state 
transition from the third excited state. The main evidence for this is the presence 
of the peak at 4-3 MeV which is identified as the one-quantum escape peak of a 
4°8MeV y-ray. Comparison with Figure 7 indicates that a 4-3 MeV y-ray 
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would not produce such a full energy peak. The lack of structure in the 3-5 MeV 
region of the large crystal spectrum precludes any definite identification of a 
4-8 MeV y-ray in this case. There is some evidence in both the small and large 
crystal spectra for the presence of a 1-75 MeV y-ray which must precede the 
4-75 MeV transition. However, the spectra in this region cannot be regarded 
as very reliable since, in order to obtain pulses of this magnitude, the necessary 
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Fig. 11.—y-Ray spectrum from 6:50 MeV state (large Nal(T1) crystal). 


discriminator setting in the y-channel was low enough to admit a large number of 
background pulses from the 128I activity in the erystal. Hence the absence of a 
sharp peak at 1-75 MeV cannot be taken as definite evidence against the presence 
of a 1-75 MeV y-ray. cs 

From these coincidence y-ray spectra we conclude that the principal y-ray 
- transitions between the low-lying levels of C are those shown in Figure 12 below. 


V. DISCUSSION 
(a) Energies of Baxcited States of 40 

Our time-of-flight measurement of the energy of neutrons to the oS 50 MeV 
state of “C is our most precise estimate of the excitation energy of an "C level. 
Our value of 6:53 -+0-02 MeV is in agreement with the values of 6-50 +0-03 cat 
found by Bent e¢ al. (1955), and 6-52 +005 MeV by Sample et pee ee 
y-ray energy measurements, but is not in agreement with i e value o 
6-40-+0-04 MeV found by both Johnson (1952) and Cerineo (1956) eerie fie 
ments of the energy of this same neutron group, using photographic : ‘ st nor 
with the value of 6-476-L0:020 MeV found by Neilson, Dawson, an. : 0 hae 
(1959), also using a time-of-flight technique, although the disagreement, in thi 


case, is not large. The accuracy of our determination depends only on the 
; ; 
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accuracy of measurement of the neutron flight path and of calibration of the 
variable delay. Flight paths could be readily measured to better than 0-5 per 
cent. and the variable delay was known at all settings to within 0-5 mys, by 
comparison with the period of a BC221 oscillator. We have been unable to 
detect any systematic error in our measurements to account for these discrepancies 
between our result and those of other workers. 


The other estimates of excitation energies from our results are the value of 
4-2+0-1 MeV for the 4-26 MeV state, from a y-ray measurement, and an energy 
difference of 0:50-+0-05 MeV between the 4-26 and 4-75 MeV states from a 
delay distribution of neutrons. Both of these values are in agreement with those 
given by Ajzenberg-Selove and Lauritsen (1959), but are not sufficiently precise 
to influence their values in any way. 


(b) y-Rays 

Since the transition probabilities for the y-rays observed are not among those 
calculated by Kurath (1957) on the intermediate coupling model, we have been 
tempted to apply the cruder estimates of Weisskopf (1951) to find the most 
probable multipolarities of the observed transitions. In doing this we have been 
encouraged by the work of Lane and Radicati (1954) on A=13 nuclei, which 
showed that the intermediate coupling and experimental values of the transition 
probabilities differed from the Weisskopf estimates by a factor of only about 10. 
Table 3 lists the Weisskopf transition probabilities of interest. 


TABLE 3 
WEISSKOPF ESTIMATES FOR TRANSITION PROBABILITIES IN UC 


H T (M1) T (£2) T (M3) 
(MeV) (sec) (sec) (sec) 
0:5 opaOl 108 4 
1-75 2x 1014 6 x 101° 3 x 104 
2-25 4x 104 aoe 2x105 
2-74 exons 6 x 104 6 x 105 
4-26 3x 10% oc 10” 107 
4-75 4x 10% Oc 1032 3x10? 
6-50 Lor 5 X 1038 3x 108 


—_—» 


Since the J,=1 stripping patterns fix the parities of the ground and first four 
excited states of 40 as odd, the only possible multipole transitions to be observed 
between these levels are M1, H2, and M3. EH4 and higher order multipoles 
have mean lives too long to be observed by our technique. 


The comparable intensities of the 6-5 MeV and 1-75 or 2-25 MeV y-Tays 
in the decay of the 6-50 MeV state is readily accountable if, and only if, the 
6-5 MeV transition is H2 and the lower energy transition M1, particularly since 
at least some #2 enhancement is to be expected. Now in Section I it was pointed 
out that the most likely value of the ground state spin is 3 /2; this value will be 
assumed throughout the following discussion. In this case the spin of the 
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6-50 MeV state would be 7 /2~ which is the assignment given by Ferguson et al. 
(1958) to the corresponding level in ™B. The possibility that the 6-5 MeV 
y-Tay is M1 or M3 is very remote in view of the comparable intensity of the lower 
energy component. 

The existence of the 2-3 and 4-8 MeV y-Tays in this decay scheme would 
indicate that both the 1-75 and 2-25 Mev transitions occur, and, since both 
would have to be M1 to compete with the 6-5 MeV H2 transition, the spins of 
the 4-75 and 4-26 MeV states would be restricted to 5/2, 7/2, or 9/2. Also, if 
both the 4-75 and 4-26 MeV ground state transitions are observed they cannot 
be M3, since this would correspond to a mean life of about 50 mys which would 


(/2>53/27) 


"Cc 


Fig. 12.—Spin assignments for low excited states of "C. 


have shown up as a considerable distortion of the delay peaks. This would 
then eliminate 9/2 as a possible spin for either the 4:75 or 4°26 MeV levels. 
While the presence of the 4-26 MeV y-ray is certain, so chat He vi case these 
conclusions are valid, the evidence for the 4-7 5 MeV y-ray is certainly less ae 
However, it should be noted that the absence of this y-ray would not affect the 
argument as to the spin of any level other than that at 4-75 MeV. 

The absence of any evidence in the decay scheme of the oe 50 MeV level for 
a cascade through the 2-01 MeV level suggests that this eared ae | 
is not M1, and consequently that the spin of the 2-01 Mea level is a i 4 / 
or 9/2. However, the expected intensity for H2 2 M3 is so low ae a oF 
these would be observed, and so spin values of 1/2 or a8 He ee inca i: 
spin of this level is 3/2 then the 4:26 MeV and possible 2 se e : an : 
the decay of the 4-26 MeV state would be both M1 or both #2. Since an upp 
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limit of 24 per cent. is placed on the intensity of the 2-25 MeV component 
relative to the 4:26 MeV component, and the theoretical relative intensities for 
M1 and #2 transitions are 1/10 and 1/20 respectively, it is unlikely that both 
are M1, but both H2 is possible; i.e. if the 2-01 MeV state is 3/2-, then the 
4:26 MeV state would be 7/2-. However, if the 2-01 MeV state is 1/2- then the 
4:26 MeV state could be 5/2- or 7/2-, the unobserved 2-25 MeV transition being 
of higher multipole order than the 4-26 MeV transition in either case. Thus a 
spin of 7/2- for the 4-26 MeV state may be associated with 1/2- or 3/2~ for the 
2-01 MeV state, while 5/2- for the 4-26 MeV state is associated with only 1/2- 
for the 2:01 MeV state. All these spin conclusions are summarized in Figure 12, 
in which are also shown the y-ray transitions observed (firm lines) and inferred 
(dotted line). 

In conclusion, it is interesting to note the extent to which the spin values 
found here for C agree with those found for “B and those predicted by the 
intermediate coupling model. The assignment of 1/2- for the first excited state 
in 1B seems fairly definite, while 5/2— and 7/2- for the second and fourth excited 
states respectively are found by Ferguson ef al. (1958). Our results are in 
complete agreement with these assignments, but in the third excited state there 
is an anomaly. Ferguson et al. (1958) have observed a 12 per cent. decay of 
the third excited state through the first, from which they conclude spin 3/2- 
for the third state. However, Wilkinson and Alburger (1959) are unable to 
reconcile a spin value of less than 5/2- for this level with their measurements on 
the 6 decay of Be. Our results for “C make a spin of 3/2- for the third excited 
state of “C seem unlikely also. The intermediate coupling model requires 
spins of 5/2- and 7/2- for the second and third excited states (or vice versa), 
and whilst Wilkinson and Alburger’s results for 4B and also our results for UC 
are in accord with this requirement, it has not yet been shown whether the model 
is capable of predicting the H2 enhancement necessary to account for the cascade 
through the first excited state observed by Ferguson et al. 
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ALPHA-PARTICLES FROM THE REACTION ‘Li(p,;)*Be*(«)*He 
By D. 8. GEMMELLT 
[Manuscript received February 5, 1960] 


Summary 


The «-particle spectrum from the reaction *Li(p,y)*Be*(«)*He has been determined 
in coincidence with y-rays. The results indicate that, apart from the broad 2-9 MeV 
level, there are no additional levels between the ground state and 7 MeV with intensities 
greater than 3 per cent. of the 2-9 MeV level. It is shown that the shape of the 2-9 MeV 
level can be fitted satisfactorily by a dispersion formula if a large value for the reduced 
width of the level is assumed. 


I. INTRODUCTION 

Several states in ®Be at excitations up to 15 MeV have been reported by 
various observers ; but, so far, the only ones which seem to be present with any 
degree of certainty (Ajzenberg and Lauritsen 1959) are the narrow 0*+ ground 
state ((~5 eV), the broad 2+ level at 2-9 MeV (I'~1-2 MeV), and a broad 4+ 
level at 11:6 MeV ([('~6-7 MeV). These levels correspond well with the only 
three states expected in this region of excitation on the basis of either a simple 
shell model or an «-particle model of Be. 


Nevertheless, the evidence for other even states in ®Be at 4-1, 5-3, and 
7-5 MeV is considerable (Titterton 1954; Ajzenberg and Lauritsen 1955), and 
it was with a view to obtaining further evidence on the possible existence of 
these levels that the present experiment was performed. 


IT. APPARATUS 

Alpha-particles from the 441 keV resonance in the "Li(p,y)’Be*(«)*He 
reaction were detected in coincidence with y-rays using the apparatus shown 
in Figure 1. A 20 keV thick target of separated 7Li on a water-cooled copper 
backing was bombarded by a collimated beam of 450 keV protons from the 
Canberra 600 keV Cockcroft-Walton accelerator. The a-particles were detected 
at 90° to the proton beam by a 14 by 4in. CsI crystal which had been thinned 
down to 0-005 in. in order to avoid an excessive background count arising from 
the passage of high energy electrons. To reach this detector, «-particles from 
the target had to pass through a slit 1 mm wide, after which they were deflected 
by a magnetic field of 5400 G produced by a set of six permanent magnets of the 
magnetron type. The purpose of the magnetic field was to prevent the intense 
8-3 MeV «-particle group from the 7Li(p,«)*He reaction, and also direct light 
from the target, from reaching the detector. Furthermore, scattered protons 
were sufficiently deflected not to reach the crystal. To protect the phototube 
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from the stray field of these magnets it was necessary to introduce a 7 in. long 
by 1}in. diameter “ Perspex ” light pipe and to surround the cathode end of 
the phototube with a soft-iron shield. 

The y-ray detector used was a 5 in. diameter by 4 in. long NalI(T1) crystal 
coupled to an EMI 6099 phototube. This crystal was placed opposite the 
a-particle counter and adjacent to the “ Perspex ” window shown in Figure 1. 


SLIT ADJUSTMENT 


MAG. SHIELD 


CRYSTAL 


| en 


WINDOW 


MAG. FIELD 
LIGHT PIPE 


Fig. 1—Diagram of apparatus. 


III, EXPERIMENTAL 

Pulses from the two detectors were fed through amplifiers into a Garwin-type 
coincidence unit (Garwin 1950) with a resolving time of 0-1 ysec, which responded 
to y-ray pulses of greater than 4 MeV, and «-particle pulses of greater than 4 MeV. 
The output from this unit was used to operate the gate circuit in a 70 channel 
Sunvic kicksorter, thus enabling the coincident «-particle spectrum to be deter- 
mined. To avoid the effects of deterioration in the target after prolonged 
bombardment the target backing was moved slightly every 2 hr so that the 
beam struck an unused section. The procedure adopted was to run for an hour 
and then obtain a background count under the same conditions by inserting a 
2-4 usec delay in the y-ray channel. Then the target was raised and, before the 
next run, the gain of the «-counter checked by moving the slit inwards a short 
distance and determining the pulse height of the 8-8 MeV «-particles from the 
7Li(p,«x)4He reaction. In all, seven separate targets were used and the final 
result comprises spectra recorded over a total bombardment time of 82 hr using 
a beam current of 5 uA. Except at very low a-particle energies the background 
count was small compared to the total number of coincidences and the resultant 
spectrum obtained after subtraction of this background is oe in qisure gi 

The calibration of the «-particle energy scale was determined by substituting 
a ThO” «-source for the target and recording the spectra obtained when various 
thicknesses of aluminium foil were interposed between the source and the ae 
In this way it was also possible to determine the resolution ies Rea 
counter as a function of «-particle energy. The use of a long lig : pipe as i 
the resolution obtainable, which was found to vary between 2 per an a 
E,=1 MeV and 15 per cent. at H,=4 MeV. The detection efficiency of the 
crystal was constant between these limits of «-particle energy. 
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IV. RESULTS 
A statistical analysis of the data in Figure 2 indicates that no additional 
levels with widths less than 0-5 MeV are present in *Be at excitations between 
2 and 7 MeV with intensities of more than 3 per cent. of the total number of 
excited *Be nuclei. 


CHANNELS ———> 


1000 


COUNTS PER CHANNEL 


ALPHA: ENERGY, ———— 


Fig. 2.—Coincident «-particle spectrum with background subtracted. The full curves 
are theoretical fits to the data, assuming M1 and #2 radiation is involved. 


The 2-9 MeV state appears to have a significant high energy tail, and an 
attempt was made to account for this by fitting the data with a Wigner-Kisenbud 
single-level dispersion formula (Wigner and Hisenbud 1947) of the form 


[17-6 —BP!+11(B) 
[E+ A(L) —E P+ [1(#)/2}” 


where # is the energy of excitation in *Be and H, is a constant. The first term 
in the numerator of this expression is proportional to the probability of y-ray 
transitions from an initial state at 17-6 MeV to a final state at H MeV (l=1 for 
magnetic dipole and /=2 for electric quadrupole radiation). The variation of 


N(#)=const. 


\- 
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the functions [and A with energy cannot be neglected in such a broad resonance 
at low excitation (Thomas 1951). Both I and A were calculated from tables of 
Coulomb wave-functions (Bloch et al. 1951) assuming a value of the Be radius of 
4°48 X10- cm (Christy and Latter 1948). In order to fit the data it was 
necessary to assume in these calculations a value for the reduced width of 
y=11-9x10™ MeV cm and a value of H,=5-95 MeV. By using these values 
and by numerically folding the experimental resolution into the theoretical 
curves (a small correction in such a broad peak) the two curves shown in Figure 2 
for M1 and H2 radiation were obtained. The fit is Satisfactory and the results 
are consistent with other evidence (Boyle 1956) that mixed M1 and £2 transitions 
are involved. The value of the reduced width is considerably larger (0?=3-4) 
than the Wigner single-particle limit (3/?/2u.R) in agreement with results found 
in the B(d,a) reaction (Treacy 1953). 
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MEASUREMENTS OF CHANGES IN THE PHASE PATH OF RADIO WAVES 
REFLECTED FROM THE IONOSPHERE AT NORMAL INCIDENCE 


By R. W. E. MoNicoxt* and J. A. THOMAS* 
[Manuscript received December 4, 1959] 


Summary 


Fixed-frequency measurements were made at Brisbane, using pulse transmissions, 
of the changes in phase of radio waves received after reflection at normal incidence 
from the various ionospheric layers. Because of the homodyne detector used, such 
phase-path records have a very good signal-to-noise ratio. Also, due to the difference 
in behaviour of echoes of different polarization or from different ionospheric layers, 
the effective resolution is high. When the echoes had a well-defined phase, the rate 
of change of phase path with time was measured ; the slowest rates were for smooth 
night-time H, layers, when values as low as 1 m sec-! were sometimes found, whereas 
for the night-time F’, region the rates were usually in the range 10-40 msec—!. Records 
were also made of reflections from drifting #, clouds, and of the effect of underlying EH, 
ionization on reflections from the F,, region. 


I. INTRODUCTION 
Measurements of changes in the phase of a radio wave received after reflection 
from the ionosphere provide much information about the characteristics of the 
reflecting layer and the ionization along the path of the wave below the level of 


reflection. The changes are due to changes in the phase path P, where P is 
defined by the equation 


h 
p=2 | udh, 
0 


v. being the refractive index corresponding to the element of height dh. 


The present paper describes a simple experimental technique for recording 
such changes in phase path and gives a review of some of the results obtained at 
Brisbane (lat. 27-5 °S., geomag. lat. 35-7 °S.). Some of the applications of 
the phase-path technique mentioned in this review have already been described 
incidentally in other papers, but it was considered desirable to bring them all 


together in one paper, and to discuss them explicitly, along with several additional 
features. 


II. EXPERIMENTAL DETAILS 
The Brisbane phase-path recording equipment employed a simplification 
of the techniques used by Findlay (1951) and Jones (1953). Pulse-modulated 


transmissions were used throughout, in order to resolve as well as possible the 
various echoes received. 
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The complete phase-path recorder consisted of a normal fixed-frequency 
group-path recorder plus one extra unit, known as the “ phase-path unit bee 
block diagram of which is given in Figure 1. This unit consisted essentially of a 
“* phase-reference oscillator ” (see below) and a mixer, in which beats were 
produced between the output of the phase-reference oscillator and the echo 
pulses coming from the receiver. These beats were then amplified and differ- 
entiated, and pips of one polarity were selected and fed to the output terminal 
of the unit. An isolating stage prevented the output of the phase-reference 
oscillator from leaking back into the receiver, which could thus be used simul- 
taneously for making group-path records if desired. 


SIGNAL FROM LAST IF. 
STAGE OF RECEIVER 


GROUND — PULSE— 
GATING AND 
ISOLATING STAGE 


CONTROL WAVE 


INPUT 
GROUND PULSE h 


AND ECHOES 
SWITCHING VIDEO 
STAGE AMPLIFIER 


Fig. 1.—Block diagram of phase-path unit, used in conjunction with a fixed-frequency 
receiver, for raking phase-path records, 


TO C.R.O, 


GROUND PULSE 
ONLY 


DIFFERENTIATOR 
AND NEGATIVE — 
PULSE SELECTOR 


PHASE 
REFERENCE 
OSCILLATOR 


To record changes in phase path, it is necessary to employ a suitable reference 
with which to compare the phase of the received signal. This reference was 
derived from an independent phase-reference oscillator (P.R.O.). Following 
Jones’s technique, the P.R.O. was operated at a frequency near the intermediate 
frequency of the receiver, rather than near the actual signal frequency, since 
this allowed recording on any desired signal frequency merely by tuning the 
receiver to that frequency, no adjustments to the phase-path unit being necessary. 

For convenience in recording, it was desirable that the P.R.O. should make 
several beats with the received signal within the duration of any given echo. 
As the duration of the transmitted pulse was about 70 usec, the P.R.O. was 
operated at a frequency about 30-40 ke/s above the intermediate frequency. 
The positive frequency difference caused downward-sloping fringes on the phase- 
path record to correspond to decreasing phase path and vice versa. 

It was necessary for the phase of the oscillations generated in the P.R.O. to 
bear some fixed relation to the phase of the oscillations in the ground pulses. 
This was most conveniently done by causing the initial phase of the oscillations 
of the P.R.O. to be exactly locked to the phase of the intermediate frequency 
signal at the start of each ground pulse. The P.R.O. was switched on at the 
instant the transmitted pulse started and a phase-locking signal, derived from 
the ground pulse, was injected into the oscillator tuned circuit, the feedback in 
the oscillator circuit being adjusted so that, after each switch-on, oscillations 
started gradually. The P.R.O. was switched off after all wanted echoes had been 
received in each pulse-repetition cycle. It will be noted that, when synchroniza- 
tion of the P.B.O. is achieved in this way, it does not matter whether the trans- 
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mitter oscillator is itself pulse modulated or whether it runs continuously and 
only subsequent amplifier stages are pulsed. 

A check on whether perfect synchronization has been achieved can be made 
by inspection of the fringes within the ground pulse. Whenever these fringes 
appear on the records, they should take the form of a series of straight lines 
parallel to the time axis of the record. 

The output pips from the phase-path unit were fed to a cathode-ray display 
tube, where they produced brightness modulation. The pattern on the cathode- 
ray tube was photographed on continuously moving 35 mm film. A large enough 
rate of film advance was used to ensure that, even when the phase path of an 
echo was changing rapidly, the successive fringes on the record could be readily 
resolved when the film was examined in a viewer. A film speed of about 
40 cm/hr was found satisfactory when the signal frequency was near 2 Mc/s, 
with proportionately higher speeds for higher frequencies. 
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Fig. 2.—Phase-path record showing two , reflecting regions, at 95 and 110 km respectively, a 

double-hop £, reflection at 190 km, an F, reflection at 210 km, an (F, +4,) reflection at 305 km, 

and a double-hop F’, reflection at 420km. Note the very slow rate of change of phase of the 
95km E, reflection. 


III. OUTLINE OF RESULTS AND DISCUSSION 
(a) Phase Coherence and Phase Continuity 

(i) Description of These Phenomena.—A phase-coherent echo is one in which 
the phase remains constant throughout the duration of the echo, i.e. the phase 
difference between the echo signal and a purely sinusoidal reference signal at 
the same frequency would be constant. In the case of a phase-incoherent echo, 
on the other hand, rapid, virtually discontinuous jumps in phase occur during 
the period of reception of the echo. 

A good example of a phase-coherent echo is given in Figure 2, in the case of 
the echo from the Z, region at 95 km, whereas Figure 3 shows highly incoherent 
echoes in the case of the H, and 22, echoes at 95 and 190 km respectively. 


——<— 
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Another important property of echoes recorded with phase-path equipment 
is the phase continuity. This is determined by the faithfulness with which the 
phase at a given portion of any particular echo is repeated from one pulse- 
repetition cycle to the next (after an interval of 20 msec). The records show 
that all echoes, whether coherent or incoherent, show at least short-term phase 
continuity. Thus phase-coherent echo traces consist of well-defined fringes ; 
and even phase-incoherent echo traces show short, slightly irregular, lines, 
approximately parallel to the time axis, the individual lines lasting for something 
between a few seconds and half a minute. 


(ii) Cause of Phase Incoherence.—Phase incoherence occurs because the echo 
is the resultant of a number of randomly phased components, coming from 
scattering centres situated at varying ranges from the recording equipment and 
hence being received at varying time delays after the transmitted pulse. At the 
start of the received echo, only those components coming from the nearest 
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Fig. 3.—Phase-path record showing highly incoherent H, reflections at 95 km (with a double hop 
at 190 km),with quite coherent F’,-region reflections at 270 km and (27, —E,) reflections at 455 km. 
The incoherent trace at 365 km is an (#,+#,) reflection. 


scatterers are received, but as time goes on during the period a BOTA of the 
echo, other components from the more remote scatterers begin ie COM EONS: to 
the resultant. The time required to change from the eA in yee a given 
component is ineffective to the one in which it is effective is only. the rise time of 
the received echo pulse from this particular scatterer, a time cs only about 
10 usec. If the amplitude of any particular ay pomponeny is comparable 
with the amplitude of the resultant immediately prior to its ertiyal and if we iii 
signals are in substantially different phase, there will Ee a significant a a 

jump in phase as the amplitude of the new component rises from Ae to sai its 
full value. Thus the phase of the resultant of all the components, i.e. the i) ase 
of the total received echo signal, can change substantially ps Ua a few micro- 
seconds; this explains the rapid, virtually discontinuous jumps in phase men- 


tioned above as characteristic of phase-incoherent echoes. 


124 R. W. E. MCNICOL AND J. A. THOMAS 


(iii) Discussion of Short-term Phase Continuity.—It is interesting to consider 
the probable reason why we find traces showing short-term phase continuity, 
even in the case of phase-incoherent echoes. In the case of such echoes, the 
total signal received at any instant is the vector sum of a large number of randomly 
phased components, coming from discrete, irregularly situated, scattering centres. 
The speed with which these centres are moving along the line of sight is believed 
to be only of the order of a few metres per second (McNicol 1949). The phase 
of the individual components of the total signal cannot change substantially 
in a time less than it takes the centres to move distances of the order of a tenth 
of a wavelength of the exploring wave, and this time is at least a few seconds. 
The resultant signal consequently cannot change phase substantially in less than 
several seconds, i.e. whatever the phase of an echo is at a given instant, that 
phase will in general persist for at least a few seconds, except on the special 
occasions when the amplitude is passing through zero, at which times the phase 
can reverse almost discontinuously. 


(b) Time-rate-of-change of Phase Path 
In the case of phase-coherent echoes, and only in this case, it is possible to 
make statements about the rate of change of phase path with time, e.g. the 
E,-region echo at 95 km in Figure 2 above shows a change of only one fringe in 
several minutes, the F’,-region echo (at 210 km) a change of about 4 fringes per 
minute, the (F,+#,) echo (at 305 km) a change of about 4 fringes per minute, 
and the 2F, echo (at 420 km) a change of about 8 fringes per minute.* 


Very slow rates of change—say less than one wavelength per minute—occur 
only for night-time H,-region echoes, and then only rarely. More commonly 
the rates of change are of the order of several wavelengths per minute, occasionally 
increasing to values well in excess of this, e.g. at times of F,-region sudden height 
rises (see Section III (h) (ii) below). 


(c) Correlation between Phase-path and “ Swept-gain’’? Records 

The degree of coherence or incoherence of phase-path records depends on 
the degree of regularity or irregularity of the reflecting region. We would, 
therefore, expect a close correlation between phase-path records and swept-gain 
records, i.e. records made by reducing the receiver gain periodically at a 
logarithmic rate from a very high value to a very low one (McNicol, Webster, 
and Bowman 1956). This correlation is in fact found to occur. An example is 
given, in Figure 4, of swept-gain records taken simultaneously with the phase-path 
records of Figure 3. As mentioned above, Figure 3 shows at 95 km an echo 
from the H, region with highly incoherent phase-path fringes. On the corres- 
ponding swept-gain record the EH, echo takes the form of “ triangular ” patches 
which show a considerable spread in virtual height at the start of each period, 
where the receiver gain is high. Both these records indicate a high degree of 


* It was usually observed that multiple-hop echoes, such as the 2F, mentioned above, showed 
phase changing at approximately twice the rate of the corresponding single-hop case, and combina- 


tion echoes, such as the (F. 2+,), showed changes corresponding to the sum of those present on the 
component echoes. 
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irregularity of the reflecting layer. On the other hand, the echo from the F 
oe at 270 km in Figure 3, shows highly coherent phase-path fringes ana 
t ese correspond with the almost “ rectangular ” swept-gain patches at 270 km 
ie oes 4, Seah show almost no spread in virtual height, even at high gain 
n this case the indications from both records are of i 
ee reflections from a very 


(d) Simultaneous Phase-path and Group-path Recording 
It will be noted that the bottom edge of the recorded phase-path trace gives 
the pane of the group path for the working frequency, albeit with slightly less 
precision than a normal group-path record would. Phase-path records thus 
offer the possibility of conveniently comparing changes in group path with 
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Fig. 4.—Swept-gain record corresponding to the phase-path record of Figure 3 showing reflections 
from irregular H, region at 95 km (with a double hop at 190 km), reflections from a smooth F', 
region at 270km, and weak irregular (F,+H,) reflections at 365 km. 


changes in phase path for any particular echo. Such comparisons are particularly 
valuable in determining, for example, whether observed changes in group path or 
phase path are due to real changes in range of the reflecting region, or to changes 
in the ionization density along the path, below the point of reflection. 


(e) Improvement of Signal-to-noise Ratio 

It was often observed that the phase-path records were more successful in 
recording very weak echoes than were group-path records made simultaneously, 
even when both recorders used the same transmitting and receiving systems. 
This was due to the fact that, for the phase-path records, the mixer stage in the 
phase-path unit was supplied with a relatively large signal from the phase- 
reference oscillator as well as with the output from the final intermediato 
frequency amplifier in the receiver. The output from a phase-reference oscillator 
is phase coherent with all of the components in the frequency spectrum of the 
wanted echo, but not with the noise pulses, which are in random phase. We 
thus have virtually a homodyne detector. Goldman (1948) points out that an 
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improvement in signal-to-noise ratio is characteristic of systems in which the 
coherence standard of the detector is maintained beyond its normal operating 
range. For the group-path records a linear amplitude-modulation detector is 
used and there is no corresponding improvement. 


(f) Improvement in Resolution of Overlapping Traces 
It frequently happens that traces which would be unresolved when using 
normal group-path techniques with a given pulse width may be effectively 
resolved on phase-path records, as, for example, the cases shown in Figures 2 
and 9, because of the recognizable differences in the behaviour of the phase path 
of the various components of the composite echo. 
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Fig. 5.—Phase-path record showing different behaviour of o- and a-reflections from F, region. 


(g) Independence of Phase Path of o- and x-Rays 
We may utilize the fact that the rate of change of phase path is often 
markedly different for the two magneto-ionie components, to give increased 
resolution above the normal h’t records. In Figure 5 for example the o-ray 
shows a phase-path increase whilst the w-ray simultaneously shows a decrease in 
phase path. Such ‘ herring-bone ”’ patterns are not at all uncommon. 


(h) F-Region Phenomena 
(i) Hatra Traces—Group-path records of F,-region night-time reflections 
often show spreading in range, and at times show a main trace and a number 
of discrete extra traces. Phase-path records have been made to find the rate of 
change of phase path of these traces (McNicol, Webster, and Bowman 1956). 


(ii) Sudden Height Rises—From time to time at night, travelling 
disturbances in the ionosphere cause distortion of the F,-region ionization 
contours. If a large disturbance occurs, it shows itself on group- or phase-path 
records as a fairly abrupt rise in the virtual height of the reflection level at any 
given frequency. The o-ray is affected first, followed after a short interval 
(of a few minutes) by the w-ray. If the working frequency is well below the 
critical frequency of the region, the o- and w#-rays at quiet times are normally not 


. 
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resolved at night on h’t records, so the effect of the disturbance is to produce a 
temporary separation of the o- and u-ray traces. <A typical example of the 
phase-path behaviour during such an occurrence is shown in Figure 6. MeNicol, 
Webster, and Bowman (1956) have made measurements on a number of such 
records and point out that the change in group path is always much more marked 
than the change in phase path of either the o- or G-TAys. 
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Fig. 6.—Phase-path record showing different behaviour of o- and 2-reflections from F, region 
during an F,-region sudden height rise. (Only every fifth fringe is shown.) 


(1) E,-fegion Phenomena 
(i) Isolated Clouds of Ionization.—One of the most striking types of phase- 
path record is obtained when we have reflections from isolated patches (clouds) 
of #, ionization drifting horizontally at a fixed height of about 100 km. Findlay 
(1951) has discussed such records, an example of which is given in Figure 7. 


t (MIN) 


Fig. 7.—Phase-path record showing drifting H, cloud at 120 km. 


Measurements of the rate of change of phase path agree ee ee cone ii 
i th within the limits of error, a 
f decrease and increase of slant group pa 3 
aii can thus be made of the speed of movement of such gia 
We might hope to find out something about the ionization density in H, 
clouds by observing the effect the passage overhead of such clouds has on the 


128 R. W. E. MGNICOL AND J. A. THOMAS 


phase-path records of reflections from the F, region. In particular, if the 
ionization density in the cloud were greater than that of its surroundings, we would 
expect the phase path from a stationary P’, region to show a temporary decrease 
as the H, cloud passes overhead. In practice, however, the picture is considerably 
complicated by the fact that the F’,-region echo is usually changing in phase path 
quite independently of the presence of any observable £, region. Also, on many 
of the occasions on which we can be sure that an observed H, cloud has passed 
exactly overhead (say, by direction of arrival measurements or by noting the 
behaviour of multiple H, reflections), the F, region is blanketed, and on such 
occasions we lose the F,-region echo and so cannot observe its changes of phase 
path. Accordingly, no definite conclusion could be drawn from the records as to 
whether EH, echoes are associated with any increase in total electron content in 
the # region or not. 
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Fig. 8.—Phase-path record showing overlapping EZ, clouds at 125 km. 


Figure 8 illustrates a rather more common event—at 125 km we have the 
obvious overlapping of echoes from a series of H, clouds to form a trace which on 
group-path records would give no such definite indication of structure. 


Occasionally the leading and trailing edges of such cloud echoes extend to 
ranges beyond that of the main echo, giving rise to converging and diverging traces 
on group-path records and decreasing and increasing phase paths on P,t records 
as discussed by Thomas and Burke (1956). 


Sometimes the recorded trace presents a completely jumbled appearance, 
as shown in Figure 3 above, i.e. we have phase incoherence, as discussed in (a) 
above. This jumbled appearance is probably due to the fact that a large number 


of clouds are present simultaneously and the traces due to echoes from individual 
clouds are no longer separately discernible. 


(ii) Simultaneous Occurrence of Several Different Types of Ks Region.— 
Quite frequently phase-path records were obtained which show smooth and 
irregular #, regions to be present simultaneously. Figure 2 shows a very smooth 
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HE, region at 95 km, together with a more irregular H, region at 110km. A 
more complex example is shown in Figure 9, where the EZ, echo ig clearly divisible 
into three parts, at 95, 110, and 130 km respectively. The phase-coherent 
central portion shows up again at 220 km on the double-hop path with twice the 
rate of change of phase path of the Single-hop echo. These latter facts make it 
probable that the H, region responsible for the echo at 110 km is overhead, at a 
height close to 110 km, rather than being displaced to one side and baie at a 
slant range of 110 km but at a smaller actual height. If this is true, the smooth Ei, 


region at 110 km must be completely embedded in some qui i 
oa quite separate, irregularl 
distributed, H, ionization. : ee z 


t (MIN) 


Fig. 9.—Phase-path record showing three different H, reflecting regions, at 95, 110, and 130 km 
respectively, and an Ff’, region at 280 km, 


(iii) Distinction between Sequential and Constant-height Es Regions.—It 
can be stated generally that the type of sporadic-H echo classified as ‘‘ sequential 
EH,” is consistently more phase coherent and stable than the type designated as 
‘“ constant-height H,” (McNicol and Gipps 1951). This feature has been used 
as an aid in distinguishing between these two types of H, (Thomas 1956). 


(j) Effect on F,-region Phase-path Records of Underlying Irregular Es 
Ionization 

One interesting effect is sometimes noticed when echoes from a partially 
reflecting H, layer are observed simultaneously with F-region echoes. Normally 
one would expect incoherent H, echoes to be accompanied by incoherent F echoes 
and this is often so. However, in certain cases, such as the one illustrated in 
Figure 3 above, the H, echoes are mainly incoherent, while the /’, echoes (which 
involve waves passing twice through the H# region) are phase coherent. In 
addition it will be noted that, in this case, the M echo, ie. the 2F',-H, echo, is 
also phase coherent. These ohservations suggest that, whatever the mechanism 
by which the H, region produces reflection of radio waves at normal incidence 
may be, the region on such occasions has the peculiar property of looking rough 
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and irregular when viewed from beneath, of looking relatively smooth when 
viewed from above (as in M reflections), and of being virtually invisible, in 80 far 
as variations of phase are concerned, to radio waves passing right through it. 


(k) Phase-path Changes during Penetration of a Layer , 

Although in this paper no direct experimental method of recording phase 
path as a function of radiated frequency is described, nevertheless some penetra- 
tion information has been gained by operating on a fixed frequency and by taking 
advantage of the fact that on certain occasions the ionospheric changes are such 
as to cause the penetration frequency of one or other of the reflecting layers to 
pass through the operating frequency of the phase-path equipment. The value 
of this procedure is somewhat reduced by the fact that it may take quite a few 
minutes for the complete record to be made and changes in the height of the 
reflecting layer as a whole may occur during this time, as also may changes in the 
distribution of ionization at levels well below the reflecting layer, and either of 
these changes will give rise to changes in the phase path other than those 
specifically due to the penetration of the layer. 


1618 
t (HR) 


Fig. 10.—Phase-path record showing penetration of H, region at 1610 hr and penetration of 
normal # region at 1626 hr. 


Figure 10 illustrates two aspects of such penetrations. Near the start of 
the record, at 1610 hr, there is a sudden reflection ‘‘ jump ”’ from H, to normal E£, 
and the phase-path fringes show no consistency from the one region to the other. 
At 1626 hr, we have, on the other hand, a penetration of the normal E region, 
and here the phase-path fringes tend to carry over across the penetration to the 
F region. 

Computations of phase path as a function of frequency carried out by Lander 
(personal communication) for H-region penetrations show a smooth increase of 
phase path with frequency. Correspondingly it has been noted that the fixed- 
frequency phase-path record always shows a continuing increase whenever the 
wave penetrates a region. 


IV. CoNncLUSsION 
Phase-path records are easy to make and provide a very useful supplement 
to the information obtained about the ionosphere in other ways. Their good 


signal-to-noise properties and their potentialities for the resolution of overlapping 
echoes are especially valuable. 
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EXPERIMENTAL RELATIONS BETWEEN IONOSPHERIC TRUE 
HEIGHT, GROUP HEIGHT, AND PHASE HEIGHT 


By J. A. THomAs* and R. W. E. McNIcou* 
[Manuscript received December 4, 1959] 


Summary 


Computations of ionospheric phase height at a particular frequency have been 
made using true height analyses of h’f curves. Temporal changes in these computed 
values agree well with experimentally measured changes of phase path both in the 
case of regular diurnal variations and during an eclipse. The true height analyses are 
thus shown to be experimentally reliable, and at the same time we can allot an “ absolute ” 
value to the normal phase-path records. 


I. INTRODUCTION 
On the basis of ray optics the phase path of a radio wave vertically reflected 
from the ionosphere at a true height “‘ h” is given by 
h 
P=2 [ udh, 
’ 0 
where yu. is the refractive index (at the frequency concerned) corresponding to the 
element of height dh. The group path is given by the companion expression 
h 
ee | u'dh, 
0 


where p’ is the group refractive index defined by p’=y+fdu/df. 


Both the group path and the changes in phase path at a fixed frequency 
may be experimentally determined. Estimates of true height of reflection may 
be made by various means based on measured (P’, f) curves, and we are therefore 
in a position to compare the measured values of P’ and AP with those derived 
from estimates of true height. 


II. EXPERIMENTAL MEASUREMENTS 

Phase-path measurements have been made at Brisbane at a frequency 
of 5-8 Mc/s by the method described by McNicol and Thomas (1960). Enlarge- 
ment of the resulting films enables the phase fringes to be counted starting from 
some arbitrary origin of time, a decrease of phase path of one wavelength being 
registered as —1, and an increase as +1. By this means we obtain plots of AP 
values throughout the daylight hours as indicated in Figure 1. Because of 
interference and the necessity for changing film it is not possible to follow the 
phase throughout a complete day. 


* Physics Department, University of Queensland, Brisbane. 
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Care is taken to count fringes corresponding to o-ray reflections only. The 
estimated accuracy of such a count at the end of a day is +10 fringes; this is 
much greater than any shift due to oscillator instabilities, but still represents 
a negligible error in a count of some 6000 fringes. There is considerable variation 
from day to day in such curves ; these day-to-day variations are usually associated 
with similar day-to-day variations in group path. 

Simultaneous h’, f records have been taken throughout the day as part of 
the routine Ionospheric Prediction Service recording and these are available for 
reduction to (N, h) profiles. Values of h’ at 5-8 Mc/s may be taken from either 
the phase-path films or the h’f curves. Piggott (1959) has recently demonstrated 
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Fig. 1.—Plot of phase-path changes throughout the day (Brisbane, April 9, 1959). 
The small irregularities shown here are typical. 


the necessity for correcting the values of h’ for the large receiver delays commonly 
found in h’f recorders. For the equipment used the delay is parece 
i is gi D=9—0-4W where W is the wi 
independent of frequency and is given by Vagpets : 

(in Sees of the recorded trace and D is the correction in kilometres. These 
corrections have been applied to the data analysed here. 


III. TruE HEIGHT ANALYSIS 
Using the method developed by Schmerling (1958) 0 the fos 
values of sampling frequencies given by either Ventrice and Se ee g : Lee ) 
or Wright and Norton (1959), a manual reduction of a number of h’f curv 
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peen carried out to derive the corresponding (h,f) curves. A typical example of 
the results of such a reduction is shown in Figure 2. The oscillation of the 
height values about the general trend is due to the limited sampling method 
employed, and has been discussed earlier by Piggott (1954). The best one can 
do in such circumstances is to fit a “‘ mean ”’ curve such that the sum of the areas 
between the mean curve and the oscillatory curve is balanced. It is obvious 
that great care must be exercised in using the Schmerling analysis for day-time 
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Fig. 2.—A comparison of (h,f) data derived from machine analysis 
(Duncan), and sampling analysis (Schmerling) of (h’f) data for 1440 hr 
on April 10, 1959. 


records. The height calculations must be carried out at very close frequency 
intervals if false data are not to be obtained. This difficulty does not exist 
at night except in the neighbourhood of the F, penetration frequency. 


For 15 cases the corresponding (h, f) data derived using the machine analysis 
of Duncan (1958) have also been plotted in Figure 2. It will be noticed that 
the two sets of data are in substantial agreement up to about 4 Mc/s, but beyond 
this frequency the Schmerling values of true height are consistently higher than 
the Duncan values. This is so for all 15 cases, the difference at 6 Mc/s being 
usually about 3km.* Such differences may often be neglected, but become 


important when one computes the phase and group refractive indices as functions 
of height. 


* Since this work was carried out Duncan (personal communication) has inserted a slight 
correction into his programme ; this correction brings the two methods into very close agreement. 
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IV. CALCULATION OF REFRACTIVE INDICES 
Using the methods of Whale and Stanley (1950) and Shinn and Whale 
(1952), we may compute the values of pand yw’ for a frequency of 5-8 Mc/s for 
vertical o-ray propagation at Brisbane, neglecting collisions. Using the previously 


calculated (h, f) profiles we obtain y and v’ as a function of true height at a number 
of times throughout the day. 


V. ComMPUTED Group PatH 
By plotting uy; .. as a function of height we obtain curves such as Figure 3. 


h 
The area beneath this curve is | p'dh and should thus agree with h’=4P’ if 
0 
the previous calculations have been correctly carried out. In practice the area 


is computed up to the height at which the plasma frequency ( f,) is 5-6 Me/s 
(u'=3-98). Between this height and the reflection height the electron density 
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Fig. 3.—Plot of u’;., versus true height for Brisbane at 1024 hr 
on April 8, 1959. 


is assumed to increase linearly, and the transformation is made to the variable 
t=4/(1—f2/5-8?). As Shinn and Whale have indicated, we can now obtain a 
finite answer to the value of 


eva ts.s (=0) jj 
| u'dh—=C(he-a—he-e) | ptdt, 
Rs. 


ts.6 


since p’t tends to a limiting value ast—>0. For Brisbane conditions this reduces 
to 12-1 (h;.3—hs-¢) km, and this value must be added to the previously calculated 


value of the area up to hs.¢. = 

Comparisons between the computed group-path EES alone actually 
observed can now be made. In all cases the two agree within a few kilometres. 
For those curves analysed using the Duncan method me computed values tend 
to be systematically low (3-5 km on the average). For those analysed i A 
Schmerling method, however, the computed values show rather more scatter 
put on the average agree with the observed values. 
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VI. CoMPUTED PHASE PATH 
By plotting p,;., a8 a function of height we obtain curves such as Figure 4. 
The area beneath this curve is (; udh and should thus agree with }P. We have 
no records of P as such, but only changes in P over a period of time. The com- 
puted values of P should vary throughout the day in the same way as does the 
change in P derived from the fringe count. A comparison of the two sets of 
data for the eclipse day of April 8, 1959, is given in Figure 5. 
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Fig. 4.—Plot of u.;., versus true height for Brisbane at 1024 hr on 
April 8, 1959. 
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Fig. 5.—Comparison of measured phase-path variations with 
computed values (x) of phase path for April 8, 1959. 


We thus have, from the good agreement shown here and throughout other 
days, a direct independent experimental check on the validity of the true height 
estimations using the Schmerling and Duncan methods. It should be noted 


that, whereas in the computation of [wan the (h, f) curve near the reflection 
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point (where p’ is high) must be very accurately known, for the computation 
of | udh the base height of the # region (where yp is high) is rather more important. 
It is just this region which is sometimes in doubt because of strong day-time 
absorption, and recourse has to be made to extrapolation of the h’f curve. This 
undoubtedly explains some of the larger deviations shown in Figure 5. 

By adjusting the vertical position of the AP curve to get the best least 
Squares fit to the | udh points, we can allot an “ absolute” value to the phase-path 
results, i.e. we have a measure of the total phase path. The curve of Figure 5 
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Fig. 6.—Comparison of “ absolute ’’ phase path for eclipse and control days. 


has been so fitted and we can thus allot an absolute phase of 5810 wavelengths 
at 1200 hr. The same can, of course, be done for other days with an error of 
less than 20 wavelengths (~1 km)—-we have thus regained some of the information 
lost through inability to record right through the night. 

We can now compute, using “absolute”? phases, the average absolute 
phase path throughout a series of days and use this for comparison with unusual 
effects produced on certain days, e.g. by an eclipse. Such a comparison is shown 
in Figure 6. 

Such comparisons are useful mainly in giving an idea of the integrated effect 
of any ionization changes taking place. The eclipse day effect will be further 
discussed elsewhere. 


VII. CoNCLUSION 
We have shown that computed values of ionospheric electron density as a 
function of height are experimentally reliable, and that they enable us to allot 
an absolute value to the norma! phase path records. True height computations 


appear to be valid during an eclipse period. 
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THE LIFETIME AND CELL SIZE OF THE GRANULATION IN 
SUNSPOT UMBRAE 


By R. E. LouGHHEAD* and R. J. Bray* 
[Manuscript received December 14, 1959] 


Summary 

High quality photographs extending over a period of 2 hr 11 min have been obtained 
of the granulation in the umbra of a large sunspot. In agreement with earlier work, 
the mean cell size of the umbral granulation is found to be significantly less, and the 
lifetime considerably greater, than that of the photospheric granulation. Ten per cent. 
of the umbral granules last for over 2 hr. 

The majority of the umbral granules, like the photospheric granules, show no 
detectable change in brightness, size, or shape during the periods over which they are 
observed to persist as identifiable structures. Some information is obtained about the 
modes of formation and dissolution of individual granules. 

The similarity between the photospheric and umbral granulation is pointed out, 
and the possible role played by convective processes in sunspots is briefly discussed. 


I. INTRODUCTION 

Recent observations of the detail in sunspot umbrae (Bray and Loughhead 
1959) have shown that the umbral granulation, while resembling the photospheric 
granulation in appearance, differs in lifetime and cell size. Measurement of a 
32 min 42 sec sequence of photographs showed that some umbral granules last 
for at least this period; the cell size, although derived from limited data, was 
found to be significantly less than that of the photospheric granulation. The 
present paper gives the results of the measurement of the lifetime and cell size of 
the granulation in the umbra of a large spot showing clear granular structure, 
based on a longer sequence showing a greater number of granules. 

The mean cell size is found to be 2-3 sec of arc as compared with 2-9 sec 
for the photospheric granulation. Sixty-seven per cent. of the umbral granules 
last for more than 15 min, 27 per cent. for mere than 45 min, and 10 per cent. 
for over 2hr. Like the photospheric granules (Bray and Loughhead 1958b), 
most umbral granules showed no detectable change of brightness, size, or shape 
during the periods over which they were observed to persist as identifiable 
structures. 

II. OBSERVATIONS 

The observations of the granulation in the umbra of the spot were made 
with a 5-in. photoheliograph using techniques similar to those described previously 
(Bray and Loughhead 1959). Two consecutive sequences of photographs 
were obtained, each covering a period of nearly an hour. The umbral granules 
appear quite clearly on both sequences, taken with exposures three and five 
times greater than those required for the photosphere (cf. Plate 1). Together, 
the two sequences allowed the granules to be studied over a total period of 2 hr 
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141 min. Some of the brighter granules are also visible on a 10 min sequence of 
twofold exposures taken shortly before. . 

The film, which was taken under good seeing conditions, contains nearly 
40 umbral photographs whose quality approaches the theoretical performance 
of the telescope. Plate 1 shows the detail present in the umbra of the spot. 
Plate 1 (a) is a normal photograph of the spot and the surrounding photospheric 
granulation. Plate 1 (b) shows, on twice the scale, a photograph taken through 
a 0:22mm diaphragm at the prime focus with three times the exposure of 
Plate 1 (a). Granulation appears over a large part of the umbra though none is 
visible in the left-hand portion, which is considerably darker than the rest. 


5” 


Fig. 1.—Umbral granulation pattern at 12 2™ 368. The heavy line 
represents the apparent outline of the umbra. The lower histogram 
in Figure 2 was derived with the aid of this map. 


Some of the granules are so close together that they can barely be distinguished. 
Note the bright forked structure in the left-hand portion of the umbra.* 
Plate 1 (c) shows a photograph with five times the exposure of Plate 1 (a), also 
taken through the diaphragm. The umbral granules now appear brighter, but 
the exposure is still insufficient to reveal any detail in the left-hand portion. 
(There are often large differences in intensity between different regions of an 
umbra, several intensity minima usually being present (Bray and Loughhead 
1959)). Note the detachment of the left-hand tip of the fork to form a granule 
during the interval between Plates 1 (b) and 1 (ec) (80 min). 


III. CELL S1zE 
Direct measurements of the diameters of individual umbral granules provide 
at best only a rough guide to the true dimensions owing to their dependence on 
photographic contrast and on the combined instrumental profile of telescope and 


* Such forks or loops are characteristic features of sunspot umbrae ; with the correct exposure 
they, too, often show a granular appearance. 
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atmosphere. However, measurements of the mean cell size of the pattern 
defined as the average distance between the centres of adjacent granules a 
independent of these effects, provided individual granules are actually see 
This quantity therefore provides a reliable and convenient parameter for 
characterizing the scale of the pattern. A similar parameter is often used in 
the theory of convective motion. 
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Fig. 2.—Cell size of the umbral and photospheric granulation 

patterns: the histograms give the distributions of the distances 

between the centres of adjacent granules (corrected for fore- 

shortening), measured for groups of 28 umbral and 92 photospheric 

granules. The mean cell sizes are 2-3 and 2-9 sec of arc respectively. 


The cell size measurements were derived from a very good photograph 
taken with an exposure factor of five. A map of the 28 umbral granules visible 
on this photograph is shown in Figure 1. Some of the granules near the edge 
of the umbra cannot be distinguished on Plate 1 owing to the difficulty of recording 
objects of different brightness with high contrast on a single reproduction. 

The distribution of the 52 intergranular distances measured is given in 
Figure 2, which also shows the cell size distribution of the photospheric granula- 
tion, derived from a good-quality photograph (Bray and Loughhead 19586: 
ef. Plate 1) taken with the same instrument. It is evident that the cell size 
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of the umbral granulation is significantly less than that of the photospheric 
granulation, the mean cell sizes being 2-3 and 2-9sec of are respectively.* 
Observations of higher resolving power might be expected to yield somewhat 
more symmetrical distributions, with some separations of less than 1 see of are. 


IV. LIFETIMES 

The lifetimes of the umbral granules were determined by a method similar 
to that used in determining the lifetimes of the photospheric granules (Bray 
and Loughhead 1958b). The results of the measurements, which refer to 48 
granules in all, are as follows : 67 per cent. of the granules had lifetimes exceeding 
15 min, 38 per cent. exceeded 30 min, 27 per cent. exceeded 45 min, and 10 per 
cent. exceeded 2 hr. Three particularly bright granules are visible on the short 
sequence of twofold overexposures as well as on the two long sequences ; these 
granules lasted for at least 2} hr.tf Umbral granules are therefore considerably 
longer-lived than photospheric granules, whose lifetimes, according to the best 
available estimates, are of the order of 10 min (Macris 1953; Bray and Lough- 
head 19580). 


V. CHANGES IN UMBRAL GRANULES 
The umbral granules exhibit a diversity in brightness, size, and shape 
similar to that shown by the photospheric granules. By examining their 
appearances on each photograph in the sequences, an attempt was made to detect 
changes in the individual granules during the periods over which they were 
observed to persist as identifiable objects. The results are given in Table 1. 


TABLE 1 
CHANGES IN UMBRAL GRANULES 


Type of Change No. of 
Granules 
No change te =e 5.6 34 


: increase 
Brightness 
decrease 


: increase 
Size 


decrease 
Change of shape .. 


~~ hp 


ot 


Comparison of Table 1 with the corresponding results for the photospherie 
granulation (Bray and Loughhead 1958): cf. Table 2) suggests that rather more 
umbral than photospheric granules undergo no detectable change (71 per cent. 
as opposed to 57 per cent.). However, this may be an observational effect, since 


* The significance of this result is established by the fact that the observed difference between 
the mean cell sizes, 0-6 sec, is six times greater than the standard error of the difference of the 
means, which is 4/{(0-09)?+(0-06)?}—=0-1 sec, 0-09 and 0-06 sec being the standard errors of 
the means of the two distributions. 


} The lifetimes of the longer-lived umbral granules are comparable with those of the penumbral 
filaments (Bray and Loughhead 1958a). 


LOUGHHEAD AND BRAY 
PLATE i 


LIFETIME AND CELL SIZE oF SUNSPOT UMBRAL GRANULATION 


(a) Sunspot photographed at 2" 13™ E.A.S.T. on May 11, 1959. Heliographic coordinates : 
115° W., 11° N.; heliocentric angle, 17°. 

(6) Threefold overexposure taken at 11" 46™ through a 0-22 mm diaphragm at the prime focus. 
Granulation appears over a large part of the umbra. With this exposure the penumbra and outer 
parts of the umbra are burnt out and appear white on the photograph. Note the bright forked 


structure at the left. 


(c) Fivefold overexposure taken at 1h 6™, 
exposure, which, however, is still insufficient to reveal detail in the left-hand portion. 
detachment of the left-hand tip of the “ fork” to form an umbral granule during the interval 


The umbral granules appear brighter with this 
Note the 


between (6) and (c). 


The approximate position of the 


diaphragm relative to the spot is shown by the white circle. 
The seale of (b) and (c) is twice that of (a). 
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14 59™ ; the first two photographs are threefold OEE ATE SAS) pine aly © five y a aes : 
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eranule which becomes detached from: the upper prong olive fork structure, and also the bright 
6 granule which ultimately attaches itself to the lower prong. 
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the smaller scale of the umbral pattern renders changes more difficult to detect. 
For those granules showing change, variations in brightness and size are unccr- 
related and have no tendency to occur during any particular part of the observed 
lifetime. 

Plate 2 illustrates the general Stability of the umbral granulation pattern 
over a period of nearly 2hr. Although there are apparent differences from one 
photograph to another due to seeing, some granules can easily be followed over 
the entire sequence. The “tails” attached to some of the granules at 
1 hr 18 min 15 see have no real existence but are due to poor seeing on this 
photograph. 

Changes in the umbral granules are difficult to detect during their periods 
of formation and decay, when the granules cannot easily be identified as such. 
However, a few cases of well-defined births and deaths were recorded. Some 
granules are preceded by a diffuse patch of bright material, which is difficult to 
distinguish from a granule smeared by seeing. Others originate by detachment 
from the inner ends of penumbral filaments or from the tip of a “‘ fork ” structure 
(cf. Plate 2). Dissolution sometimes occurs by a granule fading into diffuse 
bright material, which may then become dark. Sometimes a granule loses its 
identity by coalescing with a neighbouring granule, penumbral filament, or fork 
structure (cf. Plate 2). 


VI. DISCUSSION 

The view that the photospheric granules are to be identified with convection 
cells is supported by two main lines of observational evidence: the cellular 
nature of the granulation pattern (Loughhead and Bray 1959; Schwarzschild 
and Schwarzschild 1959), and the observed correlation between brightness and 
line-of-sight velocity in the photosphere (Plaskett 1954 ; Stuart and Rush 1954). 
Although no measurements have yet been made of the velocities cf umbral 
granules, the general similarity of the umbral granulation to the photospheric 
granulation suggests that it also may prove to have a convective origin. If so, 
the influence of the spot magnetic field on convection currents in its umbra may 
well be responsible for the smaller cell size actually observed (cf. Chandrasekhar 
1952a), though it is not clear how the magnetic field would affect the lifetime. 
In addition, theory indicates that the retardation of the circulatory convection 
currents by the field would reduce—but not suppress—the transport of convective 
energy from the bottom of the hydrogen convection zone to the visible layers 
of the spot. 

Biermann (1941) suggested that the upward transport of energy by convection 
in a spot would be completely suppressed by the field, while Hoyle (1949) pointed 
out that, assuming the material was constrained to follow the lines of force as 
they fan out near the surface, the convected energy would be distributed over a 
greater area than in the absence of a field. Chandrasekhar (19526), on the other 
hand, has concluded that under solar conditions the onset of thermal instability 
in the presence of a magnetic field is characterized by oscillations of increasing 
amplitude rather than by cellular convection. Now that the existence of umbral 
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granulation is firmly established and some information about its lifetime and 
cell size is known, further work on the behaviour of convection currents in a 
magnetic field under solar conditions is clearly desirable. 


VII. ACKNOWLEDGMENTS 
The authors wish to thank Mr. D. G. Norton for assistance in securing the 
observations, and Mr. H. Gillett for processing the film and making the 
enlargements. 


VIII. REFERENCES 


Birrmann, L. (1941).—Cf. Cowling, T. G. (1953).—‘‘ The Sun.” (Ed. G. Kuiper.) p. 570. 
(Univ. Chicago Press.) 

Bray, R. J., and LoucHHEaD, R. E. (1958a).—Aust. J. Phys. 11: 185. 

Bray, R. J., and LoucHHEaD, R. E. (1958b).—Aust. J. Phys. 11: 507. 

Bray, R. J., and LougHHEaD, R. E. (1959).— Aust. J. Phys. 12: 320. 

CHANDRASEKHAR, 8S. (1952a).—Phil. Mag. 43: 518-9. 

CHANDRASEKHAR, 8S. (1952b).—Phil. Mag. 43: 529. 

Hoyvip, F. (1949).—‘‘ Some Recent Researches in Solar Physics.” p. 11. (Cambridge Univ. 
Press.) 

LoucHueEaD, R. E., and Bray, R. J. (1959).—Nature 183: 240. 

Maocris, C. (1953).—Ann. Astrophys. 16: 19. 

PuasKkett, H. H. (1954).—Mon. Not. R. Astr. Soc. 114: 251. 

ScHWARZSCHILD, M., and ScHWARZSCHILD, B. (1959).—Sci. Amer. 200: 52. 

Sruart, F. E., and Rusu, J. H. (1954).—Astrophys. J. 120: 245. 


ere ree 


THE CHARACTERISTIC ELECTRON ENERGY LOSS SPECTRA OF 
ALUMINIUM-MAGNESIUM AND ALUMINIUM-COPPER ALLOYS 
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[Manuscript received January 13, 1960] 


Summary 


The characteristic electron energy loss spectra of 1500 eV electrons scattered by 
aluminium-magnesium and aluminium-copper alloys have been measured. A plasma 
and a lowered plasma energy loss were observed in the alloys, and it was found that 
surface oxidation caused the disappearance of the lowered plasma loss and the appearance 
of a new modified lowered plasma loss. Though the compositions of the evaporated 
alloy specimens were unknown, the three energy losses were observed to vary con- 

. tinuously between the values observed for the three elements. This continuous change 
can be interpreted in terms of the changing free electron density in alloys of varying 
composition. It is concluded that measurements of the characteristic loss spectra of 
alloys can differentiate between losses due to individual and collective electron excitation. 


I. INTRODUCTION 

In a previous paper, possible reasons for the many discrepancies in the 
measurements of characteristic electron energy losses in solids by different 
workers have been discussed (Powell 1960b). New measurements in this 
laboratory have shown the source of some of these discrepancies and have enabled 
positive identification of the origin of many energy losses to be made. Before 
that work was carried out, however, it had been considered that an experiment 
which could be useful in determining the energy loss mechanism would be the 
measurement of the energy loss spectra of certain metal alloys. 

It had been proposed that the dominant energy loss in a metal was due to 
the collective excitation of the conduction electrons (Pines 1956; Nozieres 
and Pines 1958, 1959). For a free-electron gas, the energy loss is hw, =h(4rne?/m)*, 
where w, is the free-electron plasma frequency and m is the density of free 
electrons. Where the conduction electrons in a metal are weakly bound and 
the outermost core electrons strongly bound, the observed energy loss would be 
expected to be close to iw, and the half-width of the loss small compared to 
iw, In a metal where the binding energies of the valence and outermost core 
electrons are similar, the observed energy loss would be displaced appreciably 
from iw, and the half-width of the loss peak would be of the order of fw, 

It had also been suggested that many of the energy losses were due to 
interband electronic transitions. Rudberg and Slater (1936) compared the loss 
spectrum of copper observed by Rudberg (1936) with one calculated from previous 
studies of the band structure, with fair agreement. More recently, Cauchois 
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(1952), Watanabe (1954, 1956), and Leder, Mendlowitz, and Marton (1956) 
compared the energy losses observed in a number of elements with the displace- 
ments of subsidiary absorption maxima from the K and LZ edges in the X-ray 
absorption spectra. These authors found a degree of correlation between the 
loss values and the fine-structure displacements, and Leder, Mendlowitz, and 
Marton found some dependence of the energy losses on the inverse square of the 
lattice constant in elements of similar structure. Viatskin (1958) has also found 
fair agreement between some observed energy losses and values calculated using 
the interband transition hypothesis. 


Measurements of the loss spectra of binary alloys of varying composition 
might therefore differentiate between these two mechanisms for the origin of a 
particular energy loss. If the dominant energy losses in two elements were due 
to plasma excitation, one might expect that a gradual shift in the energy loss 
with changing electron density would be observed in the alloy spectra. Alter- 
natively, if the losses were due to interband transitions, it might be expected 
that the losses due to one element would gradually decrease in intensity while 
the losses characteristic of the second element gradually increased in intensity in 
alloys of increasing content of the second element, with possibly a small variation 
in each loss position caused by band structure changes. Alloys of aluminium, 
magnesium, and copper were chosen for investigation in the present work, as 
the dominant loss in each element had been identified as arising from plasma 
excitation by Noziéres and Pines (1958, 1959) and had been correlated with the 
X-ray absorption fine structure by Leder, Mendlowitz, and Marton (1956). 
Further, the energy losses in aluminium and magnesium had been previously 
found to be sharp and in good agreement with the free-electron values of hoy, 
while the prominent loss in copper was broad and considerably different from 
hw, The dominant losses in each element were well separated, occurring at 
approximately 10, 15, and 20eV in magnesium, aluminium, and copper 
respectively, so that any gradual changes in the basic loss would be easily 
detectable. 

Before the alloy experiments were performed, measurements of the loss 
spectra of each element were made. The results of such measurements have 
been discussed elsewhere (Powell and Swan 1959a, 1959b; Powell 1960b) and 
were sufficient in themselves to clearly identify the origin of the losses. Never- 
theless, the alloy results to be presented are interesting and also confirm the 
previous interpretations. Some of the results from aluminium-magnesium alloys 
have already been described elsewhere (Powell 1960a). 


The measurements of the loss spectra of aluminium and magnesium show 
that the loss spectra of each element consist entirely of combinations of two 
fundamental energy losses, 10-3 and 15-3 eV in aluminium, and 7-1 and 10-6 eV 
in magnesium. The larger loss in each metal was identified as the plasma loss 
and the smaller low-lying loss with the lowered plasma loss predicted to occur 
at iw,/+/2 by Ritchie (1957). The copper spectrum was found to consist of a 
7-2 eV lowered plasma loss, a 19-9 eV plasma loss, a 27-1 eV combination loss, 
and a loss of 4:4 eV considered to be due possibly to an intraband transition. 
Further work (Powell and Swan 1960; Powell 1960b) has shown that surface 
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oxidation of the specimens caused a rapid decrease in intensity of the lowered 
plasma loss and the appearance of a modified low-lying loss of 7-1eV in 
aluminium, 4-9eV in magnesium, and 3-7 eV in copper. 


II. EXPERIMENTAL PROCEDURE 

The characteristic loss spectra of electrons scattered through 90° by the 
targets were obtained with the 127° electrostatic electron spectrometer and 
apparatus previously described (Powell, Robins, and Swan 1958 ; Powell 1960a). 
After the target chamber had been gettered by several evaporations of either 
aluminium or magnesium, the alloy specimens were prepared by heating a helical 
tungsten filament on which strips of two elements had been placed. This 
procedure resulted in the deposition of alloys of unknown composition on the 
target surface and was adopted after attempts to prepare alloys by flash evapora- 
tion had proved unsuccessful. Though alloys of unknown composition were 
prepared by this technique, investigations of the loss spectra of each element 
had shown the necessity for unoxidized specimen surfaces. A primary electron 
energy of 1500 eV was used for ail measurements. 


III. RESULTS 

After the simultaneous evaporation of either aluminium and magnesium or of 
aluminium and copper, a characteristic loss spectrum could be obtained which 
showed a prominent (plasma-type) loss and a low-lying loss. Typical loss spectra 
for the elements and the two types of alloys are shown in Figure 1. It was 
found that in no case could the observed spectra be regarded as arising from the 
superposition of the spectra of the corresponding pair of elements. It was 
therefore concluded that alloys had in fact been formed and that there was 
effectively no separate aggregation of the elements. From the consistency of 
the results it would also seem that the alloys must have been reasonably 
homogeneous. 

The prominent loss in aluminium-magnesium alloys was observed to occur 
between 10-6 and 15-3 eV, the values of the plasma loss in magnesium and 
aluminium. Similarly, the prominent loss in aluminium-copper alloys occurred 
between 15-3 and 19-9 eV, the observed plasma losses in aluminium and copper. 
It was found that the low-lying loss changed its position concurrently with the 
prominent loss; whereas the low-lying loss increased as the plasma-type loss 
increased in aluminium-magnesium alloys, it decreased in aluminium-copper 
alloys as the plasma-type loss increased. The continuous changes in the dominant 
loss and the low-lying loss between the positions observed in the elements confirm 
the previous identifications as plasma loss and lowered plasma loss respectively. 
The variations in the lowered plasma loss are plotted in Figure 2 (b) as a function 
of the observed plasma energy loss. | 

As the alloy specimens oxidized, the lowered plasma loss decreased in 
intensity and was replaced by a modified. low-lying loss, as in the elements. The 
observed modified low-lying loss is plotted against the corresponding plasma 
loss in Figure 2 (a); it is also seen to increase with the plasma loss in the 
aluminium-magnesium alloys and to decrease in the aluminium-copper alloys. 
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Fig. 1.—Representative portions of characteristic loss spectra of aluminium- 
magnesium and aluminium-copper alloys together with the loss spectra of the 
three elements, recorded with a primary electron energy of 1500 eV. Curve (a) 
is the spectrum of copper, curves (5) and (c) are spectra of aluminium-copper 
alloys, curve (d) is the spectrum of aluminium, curve (e) is a spectrum of an 
aluminium-magnesium alloy, and curve (f) is the spectrum of magnesium. The 
spectra have been drawn so that the plasma peak has the same intensity in all 
cases. 
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It was found that the width of the plasma loss peak at half maximum 
intensity varied considerably in the alloys and the measured half-width is plotted 
against plasma energy logs in Figure 3 (b). As the half-width of the peak of 
elasticaliy scattered electrons was 1 ‘9, eV, the true breadth of the plasma peak 
increased by a maximum factor of approximately 2-5 in the aluminium- 


MODIFIED L.L. LOSS(eV) 


L.L, LOSS(eEV) 


6 
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PLASMA ENERGY LOSS(EV) 


Fig. 2.—The observed modified low-lying loss (a) and the low-lying loss (6) 

plotted as a function of the corresponding plasma energy loss for aluminium- 

magnesium and aluminium-copper alloys. Straight lines have been drawn 

between the points corresponding to the energy losses observed in the elements 

and are not intended to imply a fit of the alloy data. The probable error in 
the abscissa is +0:2 eV. 


| magnesium alloys and by a maximum factor of about 9 in the aluminium-copper 
alloys. In contrast, there is only a gradual increase in the half-width of the 
lowered plasma loss, as shown in Figure 3 (a). It should be noted that the large 
errors in the half-widths of the aluminium-copper alloy loss peaks, due to the 
difficulty in locating the continuous background of inelastically Scattered 
electrons, preclude a more accurate measurement of the half-width variations. 
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IV. DISCUSSION 
It has been found that the loss spectra of aluminium-magnesium and 
aluminium-copper alloys do not simply arise as combinations of the spectra of 
the constituent elements. This observation is similar to that of Gauthé (1959), 
who has recently investigated the loss spectra of several intermetallic compounds. 
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Fig. 3.—The width at half-maximum intensity of the low-lying loss (a) and 

the plasma loss (6) as a function of the plasma energy loss. Straight lines have 

been drawn through the points corresponding to the values observed for the 
elements in (a) and are not intended to represent a fit of the alloy data. 


Further, the alloy spectra show a plasma and a lowered plasma loss which both 
vary continuously between the values observed in the elements. These gradual 
changes in the energy losses, presumably with composition, confirm the previous 
collective interpretation of the origin of the losses. 
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It was suggested in a previous paper (Powell 1960b) that the lowered plasma 
loss could. indicate the effective free-electron density in the specimens. This 
view appears to be confirmed by the observed decrease in the lowered plasma loss 
with increasing plasma loss in the aluminium-copper alloys. As stated in the 
introduction, the position of the plasma loss may be determined by considerations 
other than the free-electron density. Where the plasma loss is displaced from 
the value estimated using the free-electron density, a large increase in breadth 
of the loss peak is to be expected. The plasma loss peak in aluminium-copper 
alloys is thus observed to increase from that found in aluminium, even though 
the lowered plasma loss and the calculated free-electron density are decreasing. 
The large increase in breadth of the plasma loss peak is therefore consistent 
with this interpretation. It is also possible that part of the observed increase in 
breadth is due to the fact that the conduction electrons are not moving in a 
uniform lattice in the alloys ; variations in the interaction between these electrons 
and the two types of ions could cause a broadening of the loss peak. 


V. CONCLUSION 

Measurements of the characteristic loss spectra of aluminium-magnesium 
and aluminium-copper alloys have shown continuous changes in the positions 
of the plasma and lowered plasma energy losses found in the elements, Though 
the alloy compositions were unknown, these variations in the energy losses 
confirm the previous identifications. It would be of some interest in any future 
investigation to endeavour to find the relation between loss position and alloy 
composition. 

The present results also show that an investigation of the characteristic 
loss spectra of alloys can indicate the origin of the energy losses. Though the 
energy losses in aluminium, magnesium, and copper had been previously identified, 
the alloy measurements differentiated between the individual particle and 
collective nature of the origin of the losses. If the loss spectrum of a particular 
element was so complicated that the origin of the losses could not be otherwise 
determined, it could be useful to alloy that element with another element whose 
spectrum was known in order to distinguish between the losses of different origin. 
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IONOSPHERIC REFRACTION IN RADIO ASTRONOMY 
I. THEORY 
By M. M. Komesarorr* 
[Manuscript received January 13, 1960] 


Summary 


Expressions are derived for the apparent displacements of cosmic radio sources 
at transit resulting from ionospheric refraction. The derivation takes account of 
horizontal electron density gradients whose effects often outweigh the refraction expected 
from a spherically symmetrical ionosphere. 


Using these results together with position measurements taken at 19-7 Mc/s, the 
“total thickness” of the ionosphere has been estimated. This quantity relates the 
total electron content of a vertical column to the maximum electron density. During 
the observing period it had a value of 355km, in good agreement with moon-echo 
results. The observations also indicate that a considerable improvement in observed 
source positions may be achieved by applying the theory. 


I. INTRODUCTION 

Ionospheric refraction may seriously limit the accuracy of wide-aperture 
radio-astronomy aerials operating at low frequencies. A number of authors, 
notably Bailey (1948), Belyaev (1955), Link (1957), and Chvojkova (1958a 
and 19580), have estimated this refraction assuming that the ionosphere is a 
spherical shell concentric with the Earth. Smith (1952), on the other hand, 
has shown that horizontal density gradients lead to appreciable displacements of 
sources at the zenith, indicating that an adequate theory must take account of 
departures from spherical symmetry. 

The present calculation shows that, for a source at transit, the change of 
declination depends on the total electron content of a column of unit area through 
the ionosphere, together with the north-south gradient of this quantity. ‘The 
change of Right Ascension depends to a first order only on the east-west gradient. 

The investigation had a twofold aim. It was originally prompted by the 
need to eliminate refraction effects from galactic records taken with a Mills Cross 
operating as a transit instrument at 19-7 Mc/s. For this ib was necessary to 
relate observed displacements of discrete sources to published ionospheric sounding 
data. This was accomplished, and also information was gained about the 
structure of the ionosphere above the F-region maximum. 

This paper deals mainly with the theoretical aspect of the problem and 
includes an estimate of errors in the necessarily approximate treatment. A brief 
discussion of some experimental results shows that the derived value of iono- 
spheric “ total electron content ” is in good agreement with estimates derived by 
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other methods, and also indicates the improvement to the astronomical data 
which may be attained using the derived corrections. 

A preliminary report on this work has been given previously (Komesaroft 
and Shain 1959). 


II. DEFINITIONS OF THE MAIN SYMBOLS USED 
The more important symbols used are defined as follows : 
N, number of electrons per cubic centimetre, 
f, frequency of observation, 
f,, evitical frequency of the F region, 
aati, 
r, distance to the Earth’s centre, | 
radius of the surface of maximum electron density, 
r,, radius of the lower ionospheric bounding surface, 
o, terrestrial latitude (radians), 
®, latitude of observer, 
L, longitude to the east of the observer (radians), 
k, angle between ray tangent and radius vector, 
Z, observed zenith angle, 
6, declination, 
a“, Right Ascension, 
Yo Ym, ‘ Semithicknesses ” of regions above and below maximum density 
(parabolic model), 
d, total effective thickness; d= %(y,;,+Ym)- 


III. Ray PATHS IN THE JONOSPHERE 
In discussing ionospheric parameters a spherical coordinate system, r, 0, D, 
is used with origin at the centre of the Earth, 9 being latitude and D longitude 
measured to eastward from the observer. Records of critical frequency taken at 
different places show that the refractive index varies not only with r but also : 
with g and LZ. Thus the optical length of a ray path between points a and b is 


b 
P= | u(r, , L)ds 
a 


ew pee ewes 


b 
=[ [1 +r292-+r2L2 cos? o]!u(r, ©, L)dr, 


where 9=do/dr and L=dL/dr. 


From Fermat’s principle, the differential of P with respect to changing ray 
paths is zero, and this leads to the following ray equations : 


d/ ,do\ dsfdou gale’ 

ay Un ap las Be — prs, Sin @ cos ol; ak eho Ne (1) 
d dL\ ds du 

2 pac Wester As 
a COs? a5 \=a Tene (2) 


From these equations the refraction produced by any model ionosphere may be 
calculated. 
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Following Ratcliffe (1951) we will assume a density distribution of electrons 
below the maximum value V or 


N=N UA US Ns ke wes eee ae a (3) 


where y is the distance from the maximum density layer and the scale factor y 
is of the order of 100 km. A corresponding expression is adopted for the distribu- 
tion above the maximum, y,, being replaced by Yn, Where, according to the 
moon-echo work of Evans (1957), y;, is several times larger than y,,, and 

Ym +Ym =0 i 1 
r, being the radius of the maximum density layer. 


For a radio wave of frequency f, well above the electron gyro frequency 
(usually the case in radio-astronomical work), the refractive index is given by 


pamelor faa NX, ate, ae eee (4) 


where it the square of the “ plasma frequency ”’, is proportional to N. Thus 


from (3) 
Cai A I ORG) SUE a bs Mi gy who mer vague (5) 


In terms of the F’-layer critical frequency f., X,, is defined by 
Xn=felf?. 


Now, assuming that the scale factors y,,, y’, are substantially constant over 
large distances, then corresponding to horizontal density gradients we have 
refractive index gradients given by 


On —l ax ( a 
CD Agi me | op Re adeh Slee (6) 
d9 2n Op \ yn)’ 
Op —1 ox ‘ES 
x =m Se Neue eae eanaeee nee 7 
OL 2u ob (1 z) a 


In the subsequent argument it will be assumed that these partial derivatives 
of X,, are constant over small angular ranges. 


IV. THe APPARENT CHANGE OF DECLINATION 
The discussion will be limited to the case of a source at transit. Let us first 
assume that there are no east-west refractive index gradients, so that the observed 
ray lies entirely in the meridian plane. Later it will be shown that the effect 
of these gradients on observed declinations is negligible. 


(a) The Two Refraction Components 
In Figure 1, ABC is part of a ray path. The ray tangent at B, in latitude 9, 
makes an angle k& with the radius vector to the centre of the Earth O. It follows 
from the figure that the apparent declination of the source as seen from B is D, 


h 
ee DAO ks © hos ae ce eee (8) 
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It also follows that 


cos k=dr/ds, Keeeriee- eee ieee: (9) 
sin ke=rdojie, =e eens ae (10) 

and 
(tank) /r=deldr 2. =: > swe ee ee (11) 


In Figure 2, QRST is a ray from a distant source to a terrestrial observer at 7. 
The upper and lower ionospheric bounding radii are (r,,+Ym) and (7 —Ym) 
respectively. The true declination 3 of the source is indicated by the direction 


P 


Fig. 1.—Section of a ray path in the meridian plane. OP is the 
equatorial plane and angle D is the apparent declination of the 
source seen from point B in latitude 9. 


of the ray at the point Q at which it intersects the upper boundary. Denoting 
the values of angles at this level by the suffix a, it follows from (8) and (11) 
that 


gt) + id 
d= +| oe = fT Aa ot) ee (12) 


Ve 
where ® is the observer’s latitude and r, the radius of the Earth. 


The straight line ZU indicates the apparent declination 3) a8 seen from 7’. 
Denoting the angle between 7'U and the radius vector at any level by ky, it may 


Se 
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be seen that the value of 8, is given by a relation formally identical with (12), 
namely, 


mY, tan k 
s=o+|™ a ED, EC (13) 
"e 
Thus the difference A(3) between the apparent and true declination is given by 
Mm +Y py = 
A(8)=8)—3 hog het | ae eats = BO < abey 
™m—Ym 


(since k=k, when r lies between r, and r,,—y,,). The equation for ky is 


SIN ka ou eer, SINK, we ee rae ne (15) 


lIONOSPHERIC 
BOUNDARIES 


Fig. 2.—Path of a ray through the ionosphere to an observer at 7’. OP is 
the equatorial plane. (Not to scale.) 


Here Z is the apparent zenith angle of the source and the suffix m denotes the 
value of a parameter at the level of maximum electron density. To caleulate k 
we must refer to the ray equation (1). Since according to our assumption there 


are no east-west gradients, this reduces to 
Cie yee 
ar( alae do" 
Integrating and combining this equation with (9), (10), and (15) yields — 


fr Ow 
ur sin k=r,, sin Fan + { sec Bee Bee ee (16) 
m 


‘m—Ym 
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Now taking into account that u=1 at the upper boundary, it follows from 
(14), (15), and (16) that 


where 
ah i Tmt, ie a] U 
A(8) »=4%, » sin] sin k PA Geral sec k— ar, cela 
( : : ie +Ym tm—Um 09 
Vy spe E —t 
A0).=| en AAD hy ee (18) 
Tm Um Z 
and 
Sin Ky, ==1,; S10 kenl (Pz, Ui) “Soe we e pony Ss ie ees re (19) 


In substituting the value of k from (16) into (17) and (18), a simplifying 
approximation may be made. Provided k nowhere exceeds about 45°, and the 
values of dy./89 are such that the second term on the right-hand side of (16) 
is ‘small’, that is, provided 


we may neglect this term in (16), giving 
Sin 2e=P Si ky tts as ae oe eee ei (21) 


The restriction imposed by inequality (20) will be discussed subsequently. It 
may be shown that the approximation which we have made results in an error 
equivalent to an underestimate of A(5),, of the order 


Ww 


{Yn +Y nmr) |7 mt A (8) 0° 


Thus the total north-south displacement of a source at transit is expressible 
as the sum of two terms given by (17) and (18). The term A(8),,, which will be 
designated the ‘‘ wedge component ’’, depends on angular gradients of refractive 
index and vanishes in the limiting case of spherical symmetry. The term A(é), 
depends to a first order only on the radial distribution and represents the refraction 
calculated for a spherically symmetrical medium. Accordingly it is called the 
‘ spherical component’. Each of these terms is evaluated below. It is later 
shown that the wedge component is commonly predominant, indicating that 
previous treatments of the problem based on a spherically symmetrical model 
lead to erroneous results. 


(b) The Wedge Component 
Equation (17) may be simplified by a minor approximation. From (19), 
Koa=Kom 


and, since the integral term has been assumed “ small”, (17) may be restated 
thus 


— k Tmt Yy a] 
A oa SEC Ko», m U 
SOT aa rasa 7,) sec haan © nthe ae NS ne (22) 


Tm —Ym 


EE ee ee 


IONOSPHERIC REFRACTION IN RADIO ASTRONOMY. I 159 


-\n expression for sec k may be derived from (4), (5), and (21). If M is defined 
by the relation 


M = [1 a, sec? kom |? = [1 aA sec? ko,,(1 —y?/y?)\-# ovis. (23) 


then it may be shown that 


ma pea 2 Ab 9 1 
sec k=u.M sec ky, [1 +(1 —r;,/r?2) M2 tan2 Kis | PaO ee (24) 

1-50 
1-40 

5 

i 
1*30 
1*20 
1-10 
1-00 

a O-1 0-2 0-3 Or4 O°5 O26 
a 


Fig. 3.—Graph of the function 
= 3 
bis te +o)/2o4] In [(1+64)/(1—o4)]—}, 


where o=(f2/f2) sec? Keon 


Now, provided ky,,545° and oS4, where o is defined by 
o=X,,, 80C* Kons 


the value of the surd in (24) is very nearly unity throughout the range of integra- 
tion. Thus from (22) and (24) we derive the approximation 


2 2 ben He oe 
A(3) p= oe Hom | Musear. 


w 


Wee Usk) rf ye 
Substituting from (5), (6), and (23) and integrating, 
(Yn +Y;,)*(G) 8eC? Kon, (Ff) 
a= = : ag Pe auc 25 
A(8),» 3(tn +Y;,)f? do ’ ( ) 
where 
eo ae (Lae) 4 
ao)=45} sae (jaa foe ses (26) 


The parameter @(c) is shown. as a function of o in Figure 3. For values of o 
_lying between 0 and 3, @(s) lies between 1 and 1:3. 


160 M. M. KOMESAROFF 


The error of approximation in going from (17) to (25) is equivalent to an 
overestimate of about {y,,+y,,)/rn}A(S),.. Since the error mentioned in 
Section III (a) was of similar magnitude but opposite sign, the total approximation 
error is generally less than {(y,,+Ym)/1m}A(8)1o» OF One part in ten. 

The foregoing results have been derived on the assumption that the vertical 
distribution of electron density has a parabolic form. However, it may be shown 
that (25) can be generalized to cover any form of vertical profile, provided that 
this is a smoothly varying function and has a reasonably sharp upper boundary. 
For the general form we introduce the “ effective thickness ” d of the ionosphere 
given by 

aw,,= | Nar, Sho cloths: 3 scene a eee (27) 


so that for a parabolic layer 
I=3(Y mn +Ym)- 
Equation (25) then becomes 


my, ted sec? k,,, O(f2) 


SP SdiDe ade ee ee 


where 7, is the inner bounding radius of the ionosphere. 


(c) The Spherical Component 
To evaluate the integral in equation (18) we require an expression for tan k. 
From (21) and (24), 
tan k=(r,,/r)M tan k,,{1 +(1 —r?,/r2) M2 tan? Kom) ?- 


If now MM is replaced by M, its mean value for the range of integration, and the 
resulting value of tan k substituted in (18), it may be shown that this gives a 
close approximation to A(8),, 


Tm+y, 
A@).=[ mim Tn san Kom 
Tm—-Um v 


2 - ae -3 
sts 8) men] aia 8)sr mea 


etrreth ie sin a aie E sin K | ees is sin a aR i sin | 
. s 


Tn “U,, UH fn Tn yy", m—Ym 
Sal Seta k ee ee 29 
Here K is given by oe 
tan K=M tan kom 
and 
as Ak : Z 
Hoe May=~ ("" May——1— Ir May 
Ym —Ym Um 0 Vln +Ym —Um 
ele 1+o3 
=e ee + 6b oe ee (30) 


The parameter M is shown as a function of o in Figure 4. 


‘ 
} 
| 
: 


IONOSPHERIG REFRACTION IN RADIO ASTRONOMY. I 161 


The error in approximation (29) is 


jee Dee tan Kom 


re 
™m—Ym 


) 


(= e\ —4 a 
x) —— 2 é 2 |e 
1 F a (1 ') M? tan kon --M F wl (1 as in) M? tan Kn| ‘ap 


By expressing the integrand as a first-order Taylor expansion in (M—WM), it 
may be shown that the ratio of this error term to A(8),, given by (29), is condor 
ably less than (y,, Si PA Lee ‘ * 


1-407 


Fig. 4.Graph of the function M = (1/264) In [(1-+62)/(1—o4)]. 


In Figure 5, A(S), is shown as a function of k,,, according to equation (29). 
This figure also includes a graph relating k,,, to the apparent zenith angle Z (from 
(15)). Figure 5 is based on a model ionosphere for which 

Ym=0°018r,,=120 km, 

Yin =0062r,,=410 km. 
For values of y,, and y;, differing slightly from the above, A(5), may be calculated 
by assuming a linear relation with the total thickness. 

In some applications a less accurate expression than (29) is adequate. 
Thus at high frequencies and for moderate values of ko,,, the following relation 
may be derived from (29) 

ad f? 
A(8),.=— or, i: SOOM LATIN en os Las galt oer eae. (31) 
This equation, due originally to Chvojkova (1958a, 1958b), is very convenient 
for rapid numerical calculation. 
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Numerical values of A(3), may be computed by graphical integration of the 
general refraction integral 
“tan @ 
jog 


U 


taken in conjunction with equation (21). In Table 1 a number of values calcu- 
lated in this way are presented, together with corresponding values from equation 
(29). For comparison, results given by Chvojkova’s equation (31) and an 
equation due to Bailey (1948), are also quoted (the last example quoted is one 


TABLE 1 
THE SPHERICAL COMPONENT 
- | —A(8), 
Ky fe Yn Ym ne |= — . 
m 

ie (km) (km) (km) : : : Graph. 
Bailey cri ss Eqn. (29) integ) 
25° 35-9’ | 0-25 330 | 100 300 Ted 11-0’ 13’ 13-1’ 
30° 0-1875 561 198 300 19-1’ 19’ | Al ins 21-4’ 
45° 0:25 495 99 300 51-9’ 51-6’ 1° 627" 1 16" 

45° 0-25 660 165 300 Taney Pahi<G?” 12-297 i AS U 
70° 42-2’ | 0-0576 375 he) 375 ES Re | eee: Be 17 2571" 1° 26-6’ 

| i} 


which Bailey considers). It can be seen that equation (29) yields an accuracy 
five to ten times better than the other equations. However, world maps of 
critical frequency show that in moderate latitudes 0(f2)/d9 often far exceeds 
fz, 8o that from (28) and (31) 

A(3,)>A(3),. 


Under these conditions, A(S) may be calculated using the less accurate form of 
A(8), without significant loss of overall accuracy. 


(d) The Total Change of Declination 


Combining equations (28) and (31), the angle A(3) by which the observed 
declination of a source at transit exceeds its true declination is given by 


2 
A(3) = d 800? kon foes a. ee tan Kon e= (32) 


~ 2(r,-+3d]2)f2 a9 
which may be written 
_ asec? ky,, 
A(3) = + 3ajapee) sl lateness hel Days oper nce anrdes Reai (33) 


(for ko,,, d, and w@(c) see equations (15) and (27) and Figure 3). 
Equation (32) has an accuracy of about one part in ten for sources up to 
about 45° from the zenith assuming that f is not less than 24f, sec ky, This is 


not a serious restriction since radiation of frequencies less than f, sec ky, cannot 
be observed. 
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The derivation of (32) also assumes that the first term on the right-hand side 
considerably outweighs the second. When this condition is not satisfied, a 
more accurate approximation to the second term, given by (29), should be used. 
When the condition is satisfied, however, it is noteworthy that the source displaces 


10 20 30 40 50 60 


fo) 
70° 300 


° rina fe = c0o0km 


Fig. 5.—Curves showing the relationship between A(d), and observed 


2 
zenith angle Z for various values of the parameters f,/f? and Waites 
according to equations (29) and (15). The dotted line illustrates the way 
in which the figure is used to determine A(8),. In the example shown 


Z=30°, ,—1,== 300 kim, Felf2=0-3, and A(3),=—23’. 
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ment is in the direction in which f, increases. By contrast, in a spherically 
symmetrical medium the displacement is always towards the zenith. 
Finally, certain restrictions have been imposed on the magnitude of the 
horizontal gradients. From inequality (20) and equation (22), it follows that 
the first term in (32) (the wedge component) should not exceed about 1/30 
radian or 2°. In the next section the magnitude of the change in Right Ascension 
is expressed in terms of east-west gradients. Provided this change also does not 
exceed a few degrees, the foregoing derivation, which has neglected east-west 
gradients, is still valid. 


V. THE APPARENT CHANGE OF RIGHT ASCENSION 
In general the refractive index gradient has an east-west component. We 
will designate the resulting change in Right Ascension by A(«). From equation 
(2), by an argument similar to that outlined in Section IV (5), 


ee) sec 5-sec 7, SeC ky, _ O(fe) 


A(a Dr, +3d]2)f2 aL 


The angle ¢, is the latitude of the point at which the ray intersects the upper 
ionospheric boundary. ‘To sufficient accuracy sec ¢, is given by 


sec o,=sec 2 —tan k,,, tan te) me 
Provided 
cos 5 


—— 00s ky, 51 
cos O i 


(34) has an accuracy of about 5 per cent. 


VI. APPLICATIONS OF THE THEORY 

The results just derived have been applied to a number of discrete source 
observations taken with the Sydney 19-7 Mc/s Mills Cross with a beamwidth of 
1-4° (Shain 1958). 

The observing site (lat. 34°S., long. 151° E.) is well placed for calculating 
fe and its north-south gradient, since there are a number of ionospheric sounding 
stations 800-900 km apart and close to the 150th meridian. The situation is 
less favourable, however, for calculating the east-west gradient. The nearest 
stations in this direction are at Christchurch (New Zealand) and Watheroo 
(Western Australia), each more than 2000km away. Furthermore, it was 
found on comparing sounding records taken in longitude 150° E. with Watheroo 
records for the same local times, that the time derivative of f? (as used by Smith 
(1952)), was not a sufficiently accurate measure of the east-west gradient. The 
method of computation finally arrived at required data from at least three 
stations. Since there were very few days on which these were all available at 
the appropriate times, much of the information had to be based on monthly 
medians. Consequently, the east-west component estimates are considerably 
less accurate than those of the north-south component. 
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Figure 6 shows declination measurements of the source 0383A, taken between 
September and November 1957. The local times of transit were between 2330 
and 0330. It is clear that there is good correlation between observed declination 
5) and the quantity R(8), given by equations (32) and (33). (Since the source is 
close to the zenith the second term on the right-hand side of (32) may be 
neglected.) From the “ least squares ” regression line also shown on this figure, 
the “true declination” is —37°23’. The declination measured at 85 Mc/s 
is —37° 23’+3’ (Sheridan 1958), so that the two results agree within the limits of 
observational accuracy. 


From the slope of the regression line, taking vr, a8 6500 km, we calculate 
that d=355 km. This is only about 10 per cent. lower than a value derived 
by combining Evans’s (1957) moon-echo results taken in the pre-dawn period at 


— 36°30’ 


85 Mc/s POSITION 


° 20 40 60 80 100 120 140 160 180-200 
R(5)((Mc/s)* RADIAN“) 


Fig. 6.—Observed declination (epoch 1950) of the source 03S3A related 
F — 2 
to ionospheric parameters according to equation (32). R(d)=w(o)Of,/do. 


Manchester with simultaneous critical frequency measurements made at Slough. 
In an earlier report on this work Komesaroff and Shain (1959) quoted a peate of 
550 km. The lower value is based on a more careful analysis of the experimental! 
results. 

The derived value of d was used to compute Right Ascension Pusan for 
the same observations, and Figure 7 (a) shows the observed and corrected Qe Ae suit: 
of 0383A. The improvement in declination is very marked — in mains 
Ascension is by no means as good, although the mean error in Right Ascension. 
has been reduced from 0™-4 to less than 0™-1. 

Figure 7 (b) is a similar diagram for the sOulee 09S1A based on peer One 
taken between the beginning of October and mid December 1957. oe this. 
period the source was observed within vay and a nal hours of Mn, ne 
ionospheric conditions were changing rapidly, making the period Sek or 
refraction corrections. Nevertheless, although the corrected declinations s ow 
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considerable “scatter”, their mean is within 5’ of the 85 Me/s position (Mills, 
Slee, and Hill 1958). On the other hand, there is still a systematic displacement 
in Right Ascension from the high frequency position. In the absence of more 
detailed ionospheric data, it is not possible to assess the significance of this shift. 


ee O9SIA 
etn | Oi © | 


O3S3/. 


— 36°30’ 50’ 


[o) 


| 
Onno: 
° 


DECLINATION 
N 
Q 
° 
° 
° 


w 
S 
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WwW 
(eo) 


oshja™ “2o™ 217 \a5mea5™ ooi6™ 177 187 197 29m 
RIGHT ASCENSION RIGHT ASCENSION 


(a) (b) 


Fig. 7.—19-7 Me/s positions (epoch 1950) of two sources : (a) 03S3A, (b) O9SIA. 

Observed positions are shown by open circles and positions obtained after 

applying corrections for ionospheric refraction are shown by full circles. The 
rectangles indicate 85 Me/s positions. 


VII. CoNCLUSIONS 

The main results of the present paper are expressed by equations (32) 
and (34), relating the displacement components of a radio source at transit to 
ionospheric parameters. In order to compute these displacements under com- 
monly prevailing ionospheric conditions, account must be taken, not only of the 
total electron content of a vertical column through the ionosphere but also of the 
horizontal gradients of this quantity. 

Application of the results to preliminary observations, taken mainly in the 
early morning, indicates that over several months the effective ionospheric 
thickness d was apparently constant at a value of about 355 km, in substantial 
agreement with the results of Evans (1957) for another time and place. 

. Combining this value of d with sounding data on the basis of equations 
(32), and (34), a considerable improvement in measured source positions has 
been effected. Even under the somewhat unfavourable conditions around sunrise 
time, it has been possible to achieve a ten to one reduction in mean declination 
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error. To make the declination correction it was necessary to derive detailed 
information relating to north-south gradients from the published data of several 
sounding stations along the eastern Australian coast. No such detailed informa- 
tion about east-west gradients is available, however, and the improvement in 
Right Ascension which could be achieved was correspondingly less marked. 


Since these observations were made, many more position measurements 
have been taken, and these, together with a more detailed account of the experi- 
mental procedure, will be given in a subsequent paper. 
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HYPERFINE STRUCTURE IN THE MICROWAVE SPECTRUM 
OF WATER 


Ul. EFFECTS OF MAGNETIC INTERACTIONS 
By D. W. POSENER* 
[Manuscript received March 10, 1960] 


Summary 


Further analysis of the hyperfine structure appearing in the observed microwave 
spectra of HDO and D,O, taking into account nuclear quadrupole and magnetic dipole 
interactions, leads to a complete determination of the electric field gradient tensor 
at the hydrogen nuclei in the water molecule, giving a value along the bond direction 
of (+1-59+0-04) x 101° e.s.u. 


Parameters describing magnetic effects in the observed spectra have also been 
obtained. 


I. INTRODUCTION 

The partially resolved microwave spectra of the 2,,.<-2,, transition of HDO 
at 10,278 Mc/s and of the 3,,<-2., line of D,O at 10,919 Mc/s have been previously 
described (Posener 1957, hereafter referred to as I) and the results for D,O 
analysed on the assumption that only nuclear quadrupole coupling effects 
contributed significantly to the observed hyperfine structure. It was there 
concluded that the quadrupole coupling constant of the deuteron in the direction 
of the OD bond was (eqQ)op= +353-+4 ke/s, and that it was not possible to 
account for the HDO spectrum if only the nuclear quadrupolar interaction were 
taken into account. 


Improved agreement between calculated and observed spectra has now been 
obtained by including effects due to magnetic interactions occurring within the 
molecules ; it will be shown that neglect of these effects in the previous analysis 
of D,O introduced a systematic error of some 10 per cent. in the estimation of 
(eqQ)on- 


II. THEORY 
(a) Calculation of Transition Frequencies 
Hyperfine structure in the rotational spectra of HDO and D,O arises from 
interaction of the magnetic dipole and nuclear quadrupole moments of the 
hydrogen and deuterium nuclei with the various electric and magnetic fields 
occurring in the rotating molecules. The generalized theory of such interactions 


has been published (Posener 1958, hereafter referred to as IL), and our notation 
will follow this work unless otherwise stated. 


* Division of Electrotechnology, C.S.I.R.O., University Grounds, Chippendale, N.S.W. 
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Choosing for HDO the J,JF,J,FK representation, in which I, is the spin 
of the deuterium nucleus (J,=1) and I, the spin of the nucleus of hydrogen 
(I,=4), the matrix elements of the hyperfine structure part of the Hamiltonian 
may be obtained directly from II, equation (58), noting that o—0, since 
the hydrogen nucleus has no quadrupolemoment. For conciseness we will replace 
of, ah), of), and of? by a, «(), c), and c@® respectively, with the 
rotational state energy dependence implied if not indicated by subscripts. Also, 
we use the Jz, x, subscript notation for rotational states rather than J,. Then, 
with J=2, the hyperfine structure energy matrix for HDO becomes as shown in 
Table 1, to the approximations discussed in II. 


More appropriate for D,O is the J,J,lJF'K representation, in which I, and I, 
are the spins of the two deuterium nuclei (I; =I,=1). Allowing for the molecular 
symmetry the matrix elements for /=2 and J=3, as given by II, equation (53), 
are aS shown in Tables 2 and 3. 


The magnetic dipole-dipole interaction terms «‘S) occurring in the energy 
matrices may be computed as constants from IT, equation (54), using equilibrium 
internuclear distances for sufficient accuracy since the zero-point vibrations will 
not significantly affect the results. We have used (Herzberg 1945, p. 489): 
r,=0-9584 x10-8 em, and /HOH=104° 27’. With the <J>>’s calculated from 
rotational constants given elsewhere (Posener 1953), the dipole-dipole 
terms become as shown in Table 4.* 


TABLE 4 
MAGNETIC DIPOLE-DIPOLE INTERACTION CONSTANTS 


HDO D,0 
aS —25-8 ke/s a$ 7-3 ke/s 
aS —25-6 ke/s aS —4-8 ke/s 


When a diagonalizing transformation R is applied to these matrices H, 
the resulting diagonal matrices RHR give the energy splittings of the rotational 
levels, and thus the frequency shifts from the “ unperturbed” line centres v, 
can be obtained. (The tilde denotes the matrix transpose. In the sequel an 


unprimed matrix will refer to the lower rotational state, a primed matrix to the 
upper state.) 


Although the hyperfine interactions result in considerable mixing of the 
“unperturbed ” states, so that #', (in HDO) and J (in D,O) are no longer good 
quantum numbers, it is convenient to retain them as labels for a description of 
the actual energy levels. 


* In I, p. 2, it was incorrectly stated that nuclei have magnetic moments 9.t4y 3 this should 
read g,uJ,. Because of this error the nuclear g-factor for hydrogen was taken as half its correct 
value of g,,= -+-5-585340 and the calculations on HDO were actually carried through with values 


of «(S) half those shown in Table 4. As discussed below (Section IV (c)) it is not now considered 
that a recalculation is warranted. 
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(b) Calculation of Relative Intensities 
The intensity of a transition a’<-a is proportional to |(a |p| a’)|?, where 
(a |p| a’) is the dipole moment matrix ; the total intensity A of a line is the sum 
of its components, which in the present case arise only from the degeneracy in 
M,. If wp stands for the dipole moment matrix in a representation in which the 
energy is not diagonal, then in the notation of the previous section 


AG eeRER OU stee Bete cats Se ay (1) 
and so 
A= > | RuR’ |? Ree Oe Ses Po (2) 
MyM, 


For the hyperfine structure in HDO and D,O, in the non-diagonal representations 
I,JF,1I,FK and I,1,IJFK respectively, we can use IT, equation (12), to factor out 
the dependence on F, and F (or f and F’) because the dipole moment and the 
nuclear spin operators commute. Then, omitting irrelevant quantum numbers, 
we get for the elements of p, 


u(HDO)=(JKM;ivid'K' Mp) (JF, | J’F\(EF | FF’), 
p(D,O) =(JKM uid’ K' My) (IF [PBT BA Re iutnbeer seme ticks (3) 
The transformation indicated in (1) is diagonal in F and ,, so its application 


to the expressions in (3) does not affect the first factors on the right of these 
equations. If the now transformed dipole moment matrix is put into (2) and 


the summation over M,, My, and the polarizations is carried out (Condon and 
Shortly 1953, equation 745), and if we then omit factors irrelevant to the 
hyperfine structure problem, the result can be written as 


A(HDO)=| RB(HDO)R’ |?, 


ERO) =) RBQD,OVR [9 ees. aes eee ee (4) 
where 

B(ADO)=[(2F +E), | IPE | FP’), 

BDO) =[(2F+1)E(F,F)]MLIF | ILI'B'). wees seer eee (5) 


The matrix elements B(HDO) and B(D,O) are shown in Tables 5 and 6 
respectively. 


TII. ANALYSIS OF THE OBSERVATIONS 

The preceding theory shows that it is possible to calculate the hyperfine 
structure frequencies and relative intensities for the two molecules in terms of 
the fairly small number of variables appearing in the Hamiltonian. However, 
because of the small magnitude of the interactions involved, the observed 
transitions lie close together and the experimental technique could not completely 
separate them. As shown in I, only a general profile of the intensity as a function 
of frequency was measured, with occasionally a prominent peak more completely 
resolved. Nevertheless, a considerable amount of information is contained in 
such measurements, so the problem became one of calculating line profiles and 
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adjusting the variables until adequate agreement with experiment was obtained. 
Most of the large amount of numerical work required to do this was carried out 
with the aid of the digital computer SILLIAC* in about 30 hours of computing ; 
about two-thirds of this time was, however, employed on program development 


and exploratory calculations. 


TABLE 5 
THE MATRIX B(HDO) For J=2<2 


20 
Fy,F’ a2 1,3/2 2,3/2 2,5/2 3,5/2 3,7/2 
(Po hi 
1,1/2 62 —g3t — 3? 
tone 8p 152 —(3/5)3 —(27/5)2 
2,3/2 —3t —(3/5)2 152 —(5/3)# —(56/15)2 
201 
2,5/2 —(27/5)3 —(5/3)2 (70/3)? —2/15 —4/33 
3,5/2 —(56/15)2 — 2/152 (640/21)? —(32/21)# 
sya —4/3? —(32/21)2 (288/7)3 
TABLE 6 
THE MATRIX B(D,O) FoR J=3<2 
313 
Rl Qe 0,3 23 2,4 2,5 
213 
1) 212 423 
0,2 1052 
250 

2,2 —si 302 793 
2,3 —3t (63/2)3 15/23 
2,4 —(3/2)8 (45/2)3 1653 


* Adolph Basser Computing Laboratory, University of Sydney. 
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At first a Lorentz profile was used to calculate the amount each transition 
contributed to the intensity profile at a given frequency ; however, the Doppler 
contribution to the broadening was of like magnitude, and later calculations 
took this into account with subsequent improvement in the agreement between 
calculated and observed profiles. A numerical integration method was used 
for computation of this combined effect (Posener 1959). 

It was initially hoped that the whole problem could be solved automatically 
on the computer ; a program was written to calculate the line profiles for a given 
set of initial variables, and a ‘“ steepest descent ’? method used to minimize the 
sum of squares of residuals between observed and calculated ordinates. As an 
initial approximation, values corresponding to the results of I were used, i.e. 
(eqQ)op =353 ke/s and all the c® taken as zero. Although it was possible to 
decrease the sum of squares, the results did not significantly improve agreement 
with the observed line profiles because the initial approximation was poor and 
the calculation became trapped in various shallow minima which abound in the 
multidimensional surface. Similar difficulties arose when other initial approxi- 
mations were used. 

It became clear that, because of the relatively poor initial approximation, 
the minimization of the sum of squares of the residuals was alone not a sufficiently 
powerful criterion for automatic solution of the problem. This was then broken 
into parts, an approach which proved successful. Firstly, a set of «@) and c® 
was obtained to fit the HDO profile ; secondly, a similar set was found for D,O ; 
finally, these were combined. 

In order to fit the HDO spectrum alone, «®)=434 ke/s, corresponding to 
(eqQ)op =353 ke/s, was taken as approximately correct, and trial profiles computed 
for ¢)= —500(50) +500 ke/s, c= —1000(50) +1000 ke/s, a range which ought 
to cover all reasonably expected values of these parameters. At this stage no 
distinction was made between parameter values for the upper and lower states, 
for which the difference should be quite small. It was found that general agree- 
ment with the observed profile occurred only in one very restricted area of this 
net, and by examining this region with smaller mesh sizes a good fit was obtained 
for cD)~+10 ke/s, c®~ +240 ke/s. 

At this point empirical refinements, such as adjustment of the line ‘‘ centre ” 
y, and of the line-width parameters Av, (Lorentz half half-width) and Ay, (Doppler 
half half-width) were carried out; subsequently it was found that use of the 
Doppler-Lorentz profile did not significantly affect the fit in so far as the fre- 
quencies of transitions were concerned, although the overall agreement was 
much better. For Av, the value 12-0 ke/s (calculated for the experimental 
conditions) gave satisfactory results, and trial variations of several ke/s in this 
value (keeping the total line width constant) did not appear to change the results 
oan now obtained something which looked like the observed yee 
ie. a good initial approximation, more refined techniques ae ae aes y 
examining the individual transitions and, where desirable, shifting them one 
at a time independently so as to synthesize a profile pereeine somewhat ee 
with the measured one, it became possible to “ estimate ” the frequencies of the 
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unresolved transitions with some accuracy. Following this, the effects of small 
changes in the parameters «@ and c® on the frequencies of these transitions 
gave a set of linear equations to which a least squares fit could be made in order 
to further improve the parameters. The whole process was repeated as necessary 
until significant improvement ceased to be obtained. 


TABLE 7 
HDO 2,.<-2., HYPERFINE TRANSITION FREQUENCIES AND RELATIVE INTENSITIES* 


eo 
FLF| 11/2 1,3/2 2,3/2 2,5/2 3,5/2 3,7/2 
FF 
= 17242 ~~ 
1,1/2 244-4 168-3 76-9 
14-58 12-17 2-42 
32042 a 80-41 fe 
1,3/2 322-3 246-2 154-8 80-3 265-1 
12-48 27-05 0-48 15-27 3-04 
= ae ss = 35442 
2,3/2 409-9 333-8 242-4 167-9 352-6 
2-10 1-11 47-19 0-00 7-93 
201 
41141 = _ — 33142 
2,5/2 410-2 318-8 244-3 429-0 329-7 
15-20 0-03 60-63 0-91 10-74 
— 13542 a = 14412 
3,5/2 225-6 134-2 59-7 244-5 145-2 
2-80 8-23 0-81 69-74 5:93 
16542 35042 = 
3,7/2 163-2 348-0 248-6 
. 10-79 5-88 100-00 


* The first line of each row shows the estimated frequency (bold type) and the second line 
shows the frequency calculated from the best HDO parameters (column (a) of Table 9). All 
frequencies are in ke/s above 10,278-000 Mc/s. The last line of each row gives the relative 
intensities (italics). 


Estimated (using this word in the sense described above) transition fre- 
quencies for HDO are shown in Table 7,* together with the frequencies calculated 
from the best parameters (listed in Table 9). The corresponding profiles are 
shown in Figure 1, and, to a larger scale, in Figure 2. The latter may be compared 
with the spectrogram of Figure 3, which shows the corresponding part of the 
HDO spectrum at higher gain and was obtained at the same time as the experi- 
mental results reported in I. 


* All error limits quoted for the present work are r.m.s. errors. 
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In the case of D,O, the initial approximation of pure quadrupolar coupling 
already gave a good fit with the positions of the major transitions approximately 
known (see I), and again corresponding to (eqQ)op=353 ke/s there results 
ay?— —661 ke/s and af? —263 ke/s. Since the c®) were expected to be quite 
small, a least squares fit was made directly to the estimated transition frequencies 
in a manner similar to that used for HDO, and subsequently improved as far as 
possible. 

Frequencies and parameters resulting from this analysis of D,O are shown in 
Tables 8 and 9, with the profiles drawn in Figure 4. 
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Fig. 1.—Line profiles in the 2,)<-2,; transition of HDO. The “ estimated ”’ profile was UOcanaved 
ee ie “ estimated ” frequencies of Table 7, while the “ calculated ” profile was obtained ysing 
the frequencies computed from final values of the parameters (column (a) of Table 9). 


IV. RESULTS AND DISCUSSION 
The major difficulty in the present analysis was to account for the unresolved 
structure in the HDO spectrum in the region of --90 ke/s from the line centre ; 


it can perhaps be regarded as fortunate that agreement could be obtained only 


for one specific set of values of ce) and c™ out of the initial set of 861 trial values. 
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The magnitude of cH) (~240 ke/s) gives rise to frequency shifts of up to 100 ke/s 
(Fig. 5), which is in the order of that observed in formaldehyde (Okaya 1956). 
Thus it is apparent that in any molecule in which some or all of the hyperfine 
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Fig. 2.—Line profiles in the 2.)<-2,, transition of HDO. This diagram shows part of Figure 1 
with a magnified intensity scale. 


SE===== 


TTT 


Fig. 3.—High-gain spectrogram of part of the HDO 2,,<-2,, transition. 


structure is due to both deuterium and hydrogen, both will contribute about 
equally. 


For D,O the magnetic frequency shifts are much smaller (Fig. 6), and no 
difficulty was experienced in fitting the spectrum. 
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The agreement of the final calculated frequencies with the estimated ones 
to within a few ke/s appears satisfactory. In fact, this agreement is rather 
better than might have been expected from consideration of the unresolved 
observations. 


TABLE 8 
D,0O 3,3<-2,)9 HYPERFINE TRANSITION FREQUENCIES AND RELATIVE INTENSITIES* 


on 
2,1 2,2 0,3 2,3 2,4 2,5 
2584.2 
2,0 257-7 
12-73 
2,1 303-9 362-1 
12-73 25-45 
aa _ 300+1 
0,2 300-9 359-2 300-6 469-6 
0-85 8-52 54-23 0-04 
~20 
= = = 604-41 
2,2 434-3 492-6 433-9 603-0 
0-97 9-67 2-38 50-63 
ae 41242 = 521-41 
2,3 470-7 412-0 581-1 522-3 
1-82 6-71 12-38 68-18 
= ae 46242 355-41 
2,4 352-0 521-1 462-2 354-2 
0-32 0-59 13-64 100-00 


| 

* The first line of each row shows the estimated frequency (bold type) and the second line 
shows the frequency calculated from the best D,O parameters (column (a) an Table 9). All 
frequencies are in ke/s above 10,919-000 Mc/s. The last line of each row gives the relative 


intensities (italics) 


The small r.m.s. errors quoted on the variables are those obtained by applying 
mechanically the usual theory of random errors (see e.g. Brunt 1931). The 
possible effects of systematic errors are vividly illustrated in the present problem ; 
in I for D,O it was assumed that cD)—0, leading to a good least squares fit 
(eqQ)op = +353 +4 ke/s (P.E.). The seriousness of the neglect of co) and a(® 
is shown by the values now obtained of (eq@)op= +315 +7 ke/s (r.m.s8. error), 
somewhat outside the previously estimated error limits. 

In view of the good fit of calculated and observed profiles it is now profitable 
to go on to a consideration of the results shown in Tables 8 and 9. In fact, the 


Fr 
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relatively large computational effort involved in this analysis was undertaken 
in the hope that significant results of importance to molecular theory could be 
obtained, since the possession by nuclei of quadrupole and magnetic moments 
allows them to act as test probes by means of which we can in principle measure 
at their geometrical locations the electromagnetic effects of the molecular charge 


distribution. 
TABLE 9 
PARAMETERS 
Best Calculated 
Parameters* Parameterst 
(a) (0) 
(ke/s) (ke/s) 
HDO 
Yo 10,278,245-5+1 
Avo 21-6 
Avy 12-0 
Ay, 14-2 
a 434-4410 428-3 
af? 419-1410 424-1 
oy) 6-743 
oy.) 9-043 
ae 242-046 
oy) 252-0-46 
D,O 

vo 10,919,422-5-40-6 
Avo 18-2 
Avy 12-0 
Ay, 10-0 
a? —659-547 —660-5 
af? 261-843 261-8 
ey 42-044 
cy? 20-143 


*The parameters in this column were obtained for each 


molecule separately. 


{ The quadrupolar parameters were calculated from the 


variables shown in Table 10. 


(a) Interpretation of the Quadrupole Coupling Constants 
The «@ are defined by II, equation (54), as linear combinations of 
Xee—0Y<07V /0g?>, where the g refer to principal inertial axes of the molecule. 
Because of Laplace’s equation only two of the Yee are independent. Thus 
measurements on a single molecule do not completely define the quadrupole 
coupling tensor,* and it is usual to make some simplifying assumption such as. 


* Unless certain second-order effects are observable (Kikuchi, Hirota, and Morino 1959). 


‘ 
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that the bond direction is a principal axis of this tensor, i.e. that the cross-product 
terms x,,. (g’4g) vanish. Sometimes the stronger assumption of cylindrical 
Symmetry about the bond axis is also made (as in I) so that the spectra can be 
explained in terms of only one quadrupolar coupling constant, that along the 
bond direction. 
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Fig. 4.—Line profiles in the 3,;<-2.) transition of D,O. The “ estimated ” profile was computed 
using the ‘“‘ estimated ’”’ frequencies of Table 8, while the “ calculated” profile was obtained 
using the frequencies computed from final values of the parameters (column (a) of Table 9). 


In the present problem we can obtain the Xoo for two isotopic molecules 
whose principal inertial axes are rotated with respect to one another; the 
relationship between the two sets of y,, involves the usually neglected cross- 
product terms, and thus allows of a complete determination of the tensor. There 
is, of course, an approximation also involved here, the supposition that the 
electronic distribution in the two molecules is the same, 8o giving rise to the 
same electric field gradient at the, deuterons. This is not strictly true, eee 
of the different zero-point vibrations over which the field gradient ea 
averaged (Newell 1950). However, it is known that in the methyl halides for 
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example (Simmons and Goldstein 1952), substitution of deuterium for hydrogen 
changes the observed halogen quadrupole coupling constant by less than 1 per 
cent., while Newell (1950) has shown that the calculated difference between the 
electric field gradients in D, and in HD is similarly small. Thus it is reasonable 
to assume that the electric field gradients in HDO and in D,O will be the same 
to within about 1 per cent. 
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Fig. 5.—Calculated frequencies and relative intensities in the 2,)<-2,, transition of HDO. The 
lower diagram shows the large frequency errors that occur when the magnetic effects are neglected. 


Tf &, y, and ¢ form a set of right-handed orthogonal axes with the &-direction 
along an OD bond and the y-direction in the molecular plane and pointing towards 
the bisector of the angle HOH, then, by using standard transformations together 
with a knowledge of the molecular geometry, we can express the four y,, (two 
from each molecule) in terms of three independent variables yz, yen, and ea 
which completely specify the quadrupole coupling tensor ; we find ‘ 


aY(HDO)= = 1-4996 v2, 0-74.46 72, +0 3658 y5n5 
es 1 -3769-¢¢ —0 823372, +0 +1231 Ys 
ais(D,0) = —3 -600T yz¢ +0 606072, —3 4446 yng 
a(D,O) = 1-1746y¢¢ —1 +2343 Yen +0 + 8566 Yams 
Use of these equations with the «(@) of column (a) in Table 9 then gives the planar 
components of the quadrupole coupling tensor, as listed in Table 10, and these in 
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turn give the principal values of this tensor shown in Table 11 ; back-substitution 
in the equations (6) gives the values of the «(@) shown in column (b) of Table 9. 

From the accepted value of the quadrupole moment of the deuteron 
Qp = (2-738 10-014) x 10-27 em2 (Kolsky et al. 1952; Quinn e¢ al. 1958), the 
principal values of the electric field gradient at the position of the hydrogen 
nuclei can be derived. Their magnitudes are also given in Table 11, where 
comparison is made with quantities computed theoretically by R. Bersohn 
(1960) from the best available estimate of the electron distribution in the water 
molecule (Ellison and Shull 1955). 
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Fig. 6.—Calculated frequencies and relative intensities in the 3,3<-299 transition of D,O. The 
lower diagram shows that when magnetic effects are neglected small but significant errors occur in 
calculating the frequencies. 


It will be observed that, within the experimental error, the bond direction 
is in fact a principal axis of the quadrupole coupling tensor in this molecule. 

A determination of the electric field gradient in a molecule is usually used 
to provide some information on the relation of the chemical bond to the nucleus 
under consideration. In the present case, however, the field gradient arises 
essentially from the superposition of a number of comparatively small ener ias 
and, without detailed calculation based on better information than is aie zi 
available, it is not possible to draw any significant conclusions about i . 
bond. The magnitudes of the tensor components, however, give quanti nee 
data which must be accounted for by any complete theory of electron distribution, 
and to that end they provide a guide for the theoreticians. 
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(b) Interpretation of the Magnetic Coupling Constants 

The c® are also defined by II, equation (54), and are linear combinations 
of the M® of II, equation (3). Partly because the published expression for the 
M es is incorrect and partly because of complications in carrying out the required 
vibrational averaging, it is proposed to discuss the magnetic coupling constants 
in a later paper. 

TaBLE 10 
QUADRUPOLAR COUPLING TENSOR COMPONENTS 


Kee =eQp Vee=eQp(0?V /0E?) = 315-247-7ke/s 
Xen= ep Ven=eQp(C?V /0Een) = —8-8+8-7 ke/s 
Ann= OD Vinn=ePy(0? Ven?) =—139-3+7-0 ke/s 


(c) Comparison with Other Work 
Since the completion of the above work Thaddeus and Loubser (1959) 
have used maser techniques to resolve more completely the hyperfine structure 
of the HDO transition.* Their experimental results show that the frequencies 
of lines near 10,278-150 Mc/s and 10,278-350 Me/s, unresolved here, lie further 


TaBLE 11 
PRINCIPAL VALUES OF THE QUADRUPOLE COUPLING TENSOR AND OF THE ELECTRIC 
FIELD GRADIENT TENSOR* 


Electric Fiéld Gradient Tensor 
Quadrupole Coupling (e.s.u.) 

Tensor 

(ke/s) Measured Calculated + 
Axx= 31547 oV/ox#= 1-59+4+0-04~x 1015 322d Or 
Nyy =—14047 eV /0Y*?=—0-70+0-04 x 1015 —1-49 x 1015 
Xzz=—175+10 eV /0Z? =—0-89+0-06 x 1015 —1-78 x 1015 

| | = 0-115+0-061 0-092 


* X is rotated 1° 7’+1° 10’ from the &-direction towards the bisector of the 
Z HOH (ealeulated} —9-5°). 


} Bersohn (1960): Bersohn’s a, y, z, and Z correspond to our Z, Y, X, and & 
respectively. 


from our assigned values than we had estimated ; it is for this reason that the 
footnote in Section II (a) asserts that there is no point to carrying out a recaleu- 
lation now. However, the overall agreement is very satisfactory since Thaddeus 
and Loubser find «(@)(HDO)=411-1+3-2 ke/s, (¢qQ)op=312-5 ke/s, | y |<0-15, 
e)(HDO)=12-0+1-2 ke/s, and c)(HDO)=259 -8 +2-4 ke/s, in close agreement 


*The D,O line does not seem to be suitable for this type of experimental investigation 
(Thaddeus, personal communication). 
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with the results presented here. When Thaddeus and Loubser’s values are used 
to compute a profile for HDO, we find no better agreement with our experimental 
results than that previously obtained, indicating that the major source of the 
discrepancies probably lies in our measurement technique. 


The combination of Thaddeus and Loubser’s value of «@(HDO) with the 
«(2s obtained here for D,O leads to values for Yee and yyy of about 305 ke/s. 
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CROSS SECTIONS FOR THE INTERACTION OF 14-5 MeV 
NEUTRONS WITH MANGANESE AND COBALT* 


By E. WEIGOLDT 


Cross sections for (n,p), (n,2n), and (n,«) reactions have been extensively 
investigated for neutron energies near 14 MeV. However, there are still some 
gaps, and the present investigations were undertaken to provide neutron cross- 
section data for nuclei near the proton number 28. Although manganese and 
cobalt are both monoisotopic elements and, therefore, relatively simple to study 
by activation techniques, very few measurements have, in fact, been reported. 

The target samples in these experiments were in the form of thin circular 
disks of diameter 0-5in. sandwiched between a pair of similar copper disks. 
The sandwiches were irradiated for known periods, shorter than the half-life of 
the activity to be measured, in a steady flux of (14:5+0-5) MeV neutrons 
produced by the *H(d,n)*He reaction. 

After irradiation the activated samples were removed to a heavily shielded 
scintillation spectrometer employing a 2in. high by 12in. diameter NalI(T1) 
crystal, where y-ray counting was carried out in constant geometry. For ease 
of calibration of the system only counts in the photo-peak were used. Photo- 
peak efficiencies, relative to the efficiency for annihilation quanta, were measured 
experimentally using °*Co and ??Na to give calibration points close to the energies 
of the detected y-rays. The flux at the manganese or cobalt foil was taken to be 
the mean of the front and back copper monitor foils, whose activities were 
determined by counting the annihilation quanta emanating from 9-9 min ®Cu. 
The radioactive isotopes produced were identified by their radiations and half 
lives. 

The cross sections obtained are relative to a value of (522+20) mbarn 
for the ®Cu(n,2n)®Cu reaction (the weighted mean of several published values : 
Forbes 1952; Paul and Clarke 1953; and Yasumi 1957) or, in the case of the 
longer lived activities, relative to a measured value of (1030-95) mbarn for 
the reaction ®>Cu(n,2n)*4Cu. 

Five or more independent measurements were made for each cross section 


determined, and the results given below represent the weighted mean of these 
measurements. 


Reaction ®Cu(n,2n)*4Cu 


Foils of natural copper were irradiated and the annihilation quanta from 
12-87 hr “Cu were detected and counted after the 9-9 min ®Cu monitoring 
activity had decayed. Positrons are emitted in 19 percent. of the Cu 


* Manuscript received December 9, 1959. 
{ Department of Nuclear Physics, Australian National University, Canberra. 
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Saori (Strominger, Hollander, and Seaborg 1958). The cross section 
relative to the ®Cu(n,2n)®Cu cross section of (522 +20 mbarn f 
(1030 +95) mbarn. sera gee 
This may be compared with the values of (964-+78) mbarn reported by 
Rayburn (1959) relative to a Cu(n,2n)°Cu cross section of 500 mbarn, 
(954 +130) mbarn by Paularikas and Fink (1959) relative to a ®Cu(n,2n) cross 
section of 556 mbarn, and (1087 --170) mbarn by Paul and Clarke (1953), who 
find the *Cu cross section to be (482 +70) mbarn. 


Reaction °° Mn(n,2n)54Mn 
The reaction cross section was obtained by counting the 0-84 MeV y-rays 
emitted by 291 day *4Mn, and is (825-+185) mbarn. 


Reaction °° Mn(n,a)®?V 

By detecting the 1-43 MeV y-rays following the B-decay of 3-76 min ©V, 
the cross section was found to be (27-5) mbarn. 

Paul and Clarke (1953) measured a 3-9 min B-activity and attributed the 
resulting cross section of (52+8) mbarn to the (n,«) reaction. However, the 
(n,p) reaction leads to 3-6 min 5°Cr, a pure B-emitter with end point energy of 
2-85 MeV compared to 2-5 MeV for =V. Paul and Clarke identified their 
6-activities from half-lives given in the National Bureau of Standards Circular 499 
(1952); the 3-6min *°Cr has been reported subsequently. Their observed 
B-activity therefore should be due to both the (n,p) and (n,«) reactions. Thus 
subtraction indicates a value of approximately 25 mbarn for the (n,p) reaction. 


Reaction °°Co(n,2n)*%Co 

Radioactive °8Co is formed both in the 71 day ground state and a 9 hr 
isomeric state which decays to the ground state. The ground state decays to 
58He by electron capture or positron emission, followed by a 0-80 MeV y-ray in 
99-5 per cent. of the disintegrations (Strominger, Hollander, and Seaborg 1958). 
The Co(n,«) reaction leads to 2-58 hr **Mn, whose decay scheme includes a 
0-84 MeV y-ray. By following the 0-8 MeV activity until all the °*Mn and the 
9 hr ®80o isomer had decayed, it was found that (45 +25) per cent. of the °Co(n,2n) 
reactions went through the 9 hr isomeric state. 
Measurements of the (n,2n) cross section were made by detecting and counting 
the 0-80 MeV y-activity after the decay of the 9 hr isomer. The weighted mean 
of five independent measurements gives the value of (855 +165) mbarn for this 


cross section. 


Reaction °’Oo(n,p)??Fe 

Heath, Proctor, and Reich (1959) reported that 55-6 percent. of the 
disintegrations of 45 day ®Fe were followed by a 1-10 MeV y-ray, and 44-1 per 
cent. by a 1-29 MeV y-ray. Measurement of these y-ray activities gives an 


(n,p) cross section of (80 +23) mbarn. 


Reaction °°Co(n,«)**Mn 
The y-ray activity was corrected for the presence of the 0:8 MeV y-rays 


from the (n,2n) reaction. The resuiting (n,a) cross section was found to be 
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(29-46) mbarn relative to the ®°Cu(n,2n) cross section. This is in good agreement 
with the value of (31-43) mbarn relative to a °*Fe(n,p) cross section of 110 mbarn 
reported by Blosser, Goodman, and Handley (1958). 


It is a pleasure to acknowledge the many helpful discussions the author had 
with Professor E. W. Titterton and Dr. R. N. Glover during the course of this 
work. 
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A RELATION BETWEEN IONOSPHERIC DRIFTS AND ATMOSPHERIC 
DYNAMO CURRENT SYSTEMS* 


By L. H. HEISLERt 


Travelling ionospheric disturbances have been studied at this laboratory 
over a considerable period using fixed frequency observations on 5-8 Me/s. 
They are observed as anomalies in height on ionospheric records at three stations, 
Spaced approximately 40 km apart. Munro (1950) has used normal triangulation 
techniques to determine speeds and directions and has shown (Munro 1958) 
that their average speed is 7 km/min with median directions 30° east of north 
during the winter and 120° east of north by day, during the summer months. 
In both seasons there is a marked easterly component. However, in recent 
sunspot maximum years during summer, a predominant westerly component is 
evident on certain days. 


Martyn (1955) has intimated that east-west ionization drifts in the F, 
region are produced by a north-south electrostatic field, communicated from 
the main conducting region of the ionosphere by the highly conducting path 
along the lines of magnetic force. Moreover, he suggests that there should be a 
consequent reversal in direction of drift, depending on the position of the observing 
station relative to the centres of the atmospheric dynamo current systems. This 


* Manuscript received December 16, 1959. 


{ Radio Research Board Laboratory, Electrical Engineering Department, University of 
Sydney. 
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has recently been confirmed by Rao and Rao (1958) and Skinner, Hope, and 
Wright (1958), who have found drift directions near the magnetic equator, 
towards the west during the day and towards the east at night, in complete 
contrast to observations at higher latitudes by Briggs and Spencer (1954). 

It seemed probable, therefore, that the observed reversal in east-west 
component at Sydney (33° 52’S., 151° 11’ B.) is due to a shift in position of the 
current focus relative to the observing station. The transitional latitude of 
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Fig. 1.—Solar daily variation of the north component of the 
terrestrial field at latitudes 10° apart at the equinoxes in the sunspot- 
minimum year 1902 (after Chapman and Bartels). 


this current centre can change markedly from day to day or even from hour to 
hour, as shown by Bartels (1932) or in a very complete analysis eae Bocas 
(1936). Its relative position can be determined, approximately, from ae 
expressing the quiet-day diurnal variation of the X component of the Bart 8 
magnetic field for southern latitudes 10° apart during sunspot ee ae 
1902 (after Chapman and Bartels 1940) as shown in Figure 1. It will ; ae ne 

that above 40° the X component displays a maximuin, while at high lati es 
there is a pronounced minimum in the X component curve. The transi ve 
occurs at 40°, the approximate current focus centre. If, therefore, the diurna 
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curve of X component shows a maximum, the current focus is south of the 
observing station and vice versa. 

Figure 2 is a set of curves of variation in H, the horizontal component of 
the Earth’s magnetic field at Watheroo (30° 19’S., 115° 52’ B.), for the mag- 
netically quiet days January 15, 24, 26, 1955. Since variation in declination 
is very small, these are effectively graphs of the north component of the Earth’s 
field, X. There were numerous disturbances recorded on all days and the 
median values of the east-west component of the disturbance velocity are 4-3, 
3-3, and 4-3 km/min east respectively. The curves are consistent with a current 
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Fig. 2.—Variation in the horizontal component of the Harth’s magnetic field at Watheroo for the 
magnetically quiet days January 15, 24, 26, 1955. 
—-- January 15, —— January 24, —x— January 26. 


focus travelling to the north of Watheroo and subsequently to the north of 
Sydney, since Sydney is approximately at the same geographical latitude. 
A similar set of curves in Figure 3 for the magnetically quiet days January 5, 13, 
and 14, 1957 are consistent with a current focus travelling to the south of Sydney. 
Again there were numerous disturbances on all days, but in this case median 
component velocities were 6-3, 6-4, and 4-4 km/min west respectively. These 
observed drift directions agree with those suggested by Martyn for the same 
relative positions of current focus. 

Munro (1958) also reports a marked change of directions towards the west 
after midday during the summer months of the sunspot maximum years 
1956-1957. This again is probably due to a southerly movement of the current 
focus after 1200 hr local time, due to some form of solar control only evident 
at sunspot maximum. 
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Observations show that disturbances always have a dominant north-south 
component. The east-west component, which has a diurnal reversal (Munro 
1958), can be explained by the drift hypothesis of Martyn discussed above The 
north-south component, however, reverses with season and is AWiaion 
independent of geomagnetic control. ‘ 
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Fig. 3.—-Variation in the horizontal component of the Earth’s magnetic field at Watheroo for the 


magnetically quiet days January 5, 13, 14, 1957. 
—x— January 5, —-—— January 13, —— January 14. 


Observation of east-west drift is only possible because of perturbations in 
the medium. It seems probable, therefore, that the observed directions and 
velocities of travelling ionospheric disturbances are due to perturbations travelling 
north-south in an ionized medium which is drifting east-west. 
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COMMENT ON MULTILAYER DIELECTRIC FILTERS* 
By H. L. ARMSTRONGT 


A recent note by Gascoigne (1959) on transmission in multilayer dielectric 
filters is interesting as an example of an adoption to a given problem of a result 
from an entirely different problem, namely, the Kronig-Penney model for metal 
lattices. However, one may wonder whether, in the long run, it would not be 
more economical of time, especially for students, to notice that the multilayer 
filter, the Kronig-Penney lattice, and for that matter many other problems, 
are all essentially forms of the same thing. 


That is to say, they are all cases in which there are several regions of different 
properties superimposed. The behaviour of a function is specified in each region, 
and across the boundaries between regions there is the condition that the function 
and its gradient be continuous (or, instead of the gradient, there may be something 
analogous ; e.g. in electromagnetic problems one has the quantities H and #H). 
But this situation is analogous to the cascading of electrical four-terminal net- 
works, where voltage and current are continuous across the junction between 
adjacent four-terminals. Just as the problem of cascading networks is handled 
neatly by multiplying matrices, so can the problem involving superimposed 
layers be handled also by matrices. 

In this case one would have two by eee matrices relating the values of the 
function and of its gradient (or of H and H in the electromagnetic problem) at 
each surface of the layer. There would be one such matrix for each layer, and 
the superposition of layers would be represented by multiplying the matrices. 
It has already been shown (Armstrong 1956a, 1956b) how such a representation 
can be applied to a great variety of problems. In case the layers happen to be 
of two kinds placed alternately (and the Kronig-Penny model could be con- 
sidered as such a situation) one will have results involving powers of matrices, 
and there are some special relations available to be applied to these (Pease 1952 ; 
Armstrong 1953, 1956b; Mielenz 1959). 

In conclusion, then, the writer would like to suggest that it might be more 
advantageous for students to become familiar with the general methods of dealing 
with cascaded or iterated structures by using matrices, rather than by considering 
every such problem as a special case, for which special methods have to be 
developed. 

It is interesting that, about the same time that Gascoigne’s article appeared, 
there was another, about half-way around the world (Mielenz 1959) dealing with 
much the same ERB In this latter article matrix analysis is used, as well 
as some special ways of dealing with powers of matrices. 


* Manuscript received December 22, 1959. 
ft Queen’s University, Kingston, Ontario, Canada. 
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ON THE FLOW AND EXCHANGE DIFFUSION OF ELECTRONS AND 
NEGATIVE IONS FORMED BY MOLECULAR ATTACHMENT* 


By W. G. KAUMANT 


Formation of negative ions by attachment of electrons to molecules was 
first recognized by Thomson (1916). Huxley (1959) recently discussed the 
steady-state structure of a stream of electrons and ions drifting and diffusing in 
a gas, taking account of ionization by molecular attachment. 

In general, there is a finite probability of de-ionization by detachment of 
electrons, particularly at high ratios of applied field to pressure or when attach- 
ment results in a short-lived excited state (Loeb 1955). With detachment non- 
negligible, the process may be considered as an example of ‘‘ exchange diffusion ’’t 
(Kauman 1959) and treated by a method developed for electrodiffusion (Bak 
and Kauman 1958, 1959). In the present communication we shall show that 
this method permits separation of the rates of ionization and de-ionization and 
may be readily applied to non-steady states, alternating fields, and systems with 
a variety of boundary conditions. 

Using Huxley’s notation but putting k, (instead of fv) for the rate of 
ionization and k, for the rate of de-ionization, assuming first-order reactions, 
and neglecting, for simplicity, paired production of positive ions and electrons, 


we have 
An/ot=D,Vn—W,, -Vn—kyn+kheW, eee eens (1) 


aN /at=Dyv2N —Wy:-VN—KkyN+hyn,  «.- 20s (2) 


where m=number of electrons per unit volume, 
N=number of negative ions per unit volume, 
D,=diftusion coefficient (i=n, N), 
W ,=em,H=adrift velocity, 
e=electronic charge, 
m;=mobility, 
E=electric field (antiparallel to 2-axis). 


* Manuscript received December 23, 1959. 

+ Division of Forest Products, CS R:O., Melbourne. 

{ “‘ Exchange diffusion ”* occurs when two species of particles (Z7,, My) 
M,=M,, diffuse in an external force field. 


, reacting according to 
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The system (1), (2) may be solved by first calculating the result for an 
instantaneous point source of electrons at R=0, t=0 and then integrating with 
respect to time to convert to a continuous source of strength s (electrons per unit 
time) at the origin (cf. Crank 1956). Neglecting small transient terms and terms 
above the second order in H, the non-steady-state solution for 1°Zt<1 and 


WDGAK<1 is 


[he ofl xGD_R , 
ne qa| gan? |potekevi— TT—b E jaty! 
‘Ris ot ey 12D, +x)R Se Pre 
Bie eaueat aed in T(—4, Rk 149, eA) 
where 
B=-9D,9., B=-G9DTD., 

9, =D, =(ke_D, +h Dyl |x, 9, =9,=([k,D, +k_Dyl |x, 
9.=9,4+9(B), with D(E)=[(m, —my)eE]?k,ke/x3, 
9,=9,—-H(B), 

R*=27/9,+y7/D, +27/F., R2=2/9,+y2/9,+27/F., 
A=[kmy+kym, |eH/2D,x, A=[kym, +kym,]eH/29.x, 
x=kh,-+k,. 


In the steady state, the solution is 
n= & On? R-O TOK, 4 (ARVSQ,) 
r=0 E 


+ © b,2R-O+OK, ARV (D,+x/02)}, eee (4) 


r=0 


where K,,; are modified Bessel functions of the second kind of order r+. 
The first terms of the series in (4) are 


8 
~ Arex 


ky k, 
last exp A(2—R~/Q,) + See Daiap exp Mz—R/( Q +f]. 


The solutions for NV are similar to the above, except that in the absence of a 
source of negative ions, the sign of the term in exp (dz) is negative and k, in the 
leading term is replaced by k,, with corresponding changes in the pietteerite 
in (4). Assuming an appropriate image source (cf. Huxley 1959), we can make 
n=N=0 at z=h by subtracting from each term in (3), (4), (5) a similar term with 
2 replaced by 2h—z. 
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The electronic current to an infinite plane electrode at <=h is, to a first 
approximation, 


ino == [hy thy exp AR{1—+/(1+%/229,)}]. ...... (6) 


If W=4D,k,, ki>k,, and m,>my, (5) and (6) reduce to Huxley’s equations 
(9) and (11). 

If the applied field alternates with frequency w/27, the results are of the 
same form as above, but the series terms are more complicated and depend also 
on @, and YF) (ef. equation (3)), averaged over a cycle, is 


= Ep\ 42m’ ky —Kp)?\) 


where H,/1/2=root-mean-square field. Further, A=[(m, +my)eH, cos wt] /2Z., 
and similarly for 3. 

The general equations (1) and (2) may also be solved with finite 
beundary conditions, for instance n/(0,t)=N(0,t)=40), n(h,t)=N(h,t)=1C,, 
(C,—C,)/h=g<C,. The one-dimensional steady-state solution for k,=k,=k, 
|m,—my |<k is then of the form 


newdg(h—z) --40,+ {1 --e-Az eA 
2 22 2 
eee +eAz_eAe-h) +0[(m,, —my)/I], Bae oS) 


where A~meH/D and D,m are mean values. 


It is important to note that in the systems described, the divergence of the 
electric field is finite, and the drift. velocity of the particles therefore is not 
constant. For instance, the electric field corresponding to the system described 


by equation (8) is of the order 
Erk t(efe\O,(4e—2), 02 cect cence ee. (9) 


where H,=applied field, e=permittivity. The effect of the divergence of # may 
be taken into account by perturbation procedures. In the presence of a species 
of positive ions not taking part in the exchange diffusion, the divergent term of 
the field is, to a first approximation, {eH ,(Co—C;,)/2h(O, +3C,)}(2 —3h), where C, 
is the concentration of the positive ions. The electric current in this case is 


I~ — or) 20-607 eS (Hae vee (10) 


With an alternating field, 7 is the form I(alt.)=J)+, sin wt +I, cos 2mt, where 
I, and I, are independent of CO, but depend on k, m, and w with a cut-off frequency 
given by o=2k (Kauman and Bak 1959). 


G 
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The present method not only permits a complete description of the behaviour 
of electrons and ions formed by molecular attachment, but is also amenable to 
the treatment of other electronic exchange processes, for instance transport of 
carriers in semiconductors (ef. van Roosbroeck 1950; Shockley 1950). It may 
be extended to more than two ionic species by iteration procedures. 
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MAGNETIC ANISOTROPY OF DISPERSED POWDERS* 
By F. D. StTAcEyy 


It is well known that in most rocks the ferromagnetic fraction occurs as 
small grains dispersed in a solid, magnetically inert matrix. Recently the 
magnetic anisotropy of rocks and of chondritic meteorites has been subjected 
to detailed study by the torque-meter method, and, in an attempt to obtain a 
physical understanding of the shape and crystal alignments of grains which cause 
magnetic anisotropy in these natural bodies, a number of artificial specimens 
have been prepared. Iron, nickel, and. magnetite powders were mixed into 
solidifying media and allowed to set in a cylindrical mould in a 10 kilo-oersted 
field. The torque curves of the resulting specimens reveal a surprising result. 
It appears that single crystal magnetic grains tend to string together along lines 
of force, thus producing strongly anisotropic specimens, only when the first 
magnetocrystalline anisotropy constant of the ferromagnetic material is positive. 


The Torque-meter Hxperiment 


Torque meters are generally used to measure magnetic anisotropy, and 
Bozorth (1951) has summarized the method and the results which have been 
obtained with single crystals. For measurements on rocks the simple instrument 
shown in Figure 1 is adequate. The torsional suspension 7 is a 7-5 em length 


* Manuscript received December 30, 1959. 
+ Geophysics Department, Australian National University, Canberra. 
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of 0-45 mm diameter phosphor-bronze wire, hanging from a needle NV » which can 
be turned by hand round a 360° seale H in a horizontal plane. A brass rod R 
with an adjustable mirror M clipped to it is soldered to the lower end of the 
Suspension ; cylindrical specimens are stuck to this rod with wax and hang 
freely. Only very strongly magnetic, light specimens need to be weighted. The 
torque meter is placed so that the specimens hang symmetrically in the field gap 
of an electromagnet (Stacey 1959), which is used in this experiment with 10-8 em 


N H 


Fig. 1.—Simple torque meter for anisotropy measurements 
on rocks. The magnetic field is applied perpendicular to 
the plane of the diagram. 


plane pole faces and normally operated at 10-6 kilo-oersteds over a gap of 4-0 em. 
The mirror is used with a galvanometer lamp and scale for sensitive observations 
of the angle of the specimen. ; 

To make a measurement the needle is set to a predetermined angle @ and 
the mirror adjusted to bring the light spot into the middle of the cals, The 
field is switched on and the rotation of the specimen is compensated by tuming 
the needle through an angle AO to return the light spot to its original position. 
Then AO is the angle through which the suspension has been twisted to yee: 
exactly the magnetic torque on the specimen. This process is repeated at 0=0°, 
10°, 20°,. . ., to 180°, and sometimes to 360° as a check. Then T=«A0 iat le 0 
is the torque curve of the specimen, where « is the torsional constant “4 the 
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suspension. AQ can be measured to about 0-2°, which is somewhat less than 
lper cent. of the maximum normally usable deflection. For very weakly 
anisotropic specimens the direct observation of A@ is inaccurate and a more 
Sensitive method uses calibrated deflections of the light spot without movements 
of the needle. 

For convenience the artificial specimens are set in the torque meter so that 
their known axes of anisotropy coincide with the field direction when the scale 
reading is zero. Their torque curves can then be resolved into even harmonics 
in 0: 

21(0)=—7', 810 20 7 sin 40-7 sin 60h. a. a ecse (1) 


Small disturbing terms 7, and 7, sin (0+ 9), (o representing arbitrary phase) 
may also appear. The 7, term is due to rotational hysteresis and is generally 
observable only at low fields (a specimen at §-+10° tends to “‘ remember ” the 
magnetization it acquired at 6 and thus experience a torque independent of 0). 
The 7, term may be observed with specimens having unevenly distributed 
magnetic material, suspended asymmetrically in the field. It has been noticed 
with a few rocks, but with none of the artificial specimens. 


EHaperimental Results 

It is to be expected that spherical, polycrystalline grains of ferromagnetic 
material, immersed in a non-magnetic paste and exposed to a strong magnetic 
field, will tend to align themselves in strings along the lines of force. The 
specimen which is formed when the paste has set will have a magnetic anisotropy 
which may be simulated by a pair of individually isotropic ball-bearings stuck 
together as specimen 1 (see Table 1). The similarity of the torque curves of 
the other specimens to that obtained with specimen 1 is a strong indication that 
the observed anisotropies were principally due to the stringing together of 
magnetic particles. 

However, in spite of the similarity of shape of all the torque curves, their 
amplitudes differ most remarkably. The specimens of nickel and magnetite 
powder, which were mixed with about 1000 times their own volume of plaster of 
paris (5, 6, and 9), wetted, and allowed to set in the field for about 5 hr, exhibited 
very little or no anisotropy, whereas the iron spec'mens (2, 3, and 4), prepared 
in the same way, were strongly anisotropic. 

We may consider the behaviour of a specimen containing single (cubic) 
crystal magnetic particles which rotate in the setting field to give [100] crystal 
alignment with the field and random crystal directions in the perpendicular 
plane. Neglecting the stringing effect it would have a torque curve given by 


eae eT vs acy ie OT eee 
1=K,(5 sin 20+ 55 sin 16) +K( ig sin 20+ 35 sin 40 — 18 sin 60), 


where K, and K, are the first and second magnetocrystalline anisotropy constants 
of the material. Since for iron K,~+4 x10? ergs em- and K, is negligible 


for this purpose, it is apparent that the sign of the sin 40 term in the measured 
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torque curves is opposite to that which would result from crystal alignment. 
Further, no distinction can be made between the shapes of the torque curves 
for the specimens containing iron (2, 3; and 4) and the curve for the pair of 
ball-bearings (specimen 1), so that crystal alignment does not contribute to the 
anisotropy of the iron specimens. This is expected from the polycrystalline 
nature of carbonyl iron powder (Pfeil 1947). 


The torque exerted on specimen 3 with its axis of anisotropy at a fixed angle 
(30°) to a variable field H was measured over the range 3-9-11-5 kilo-oersteds. 
It showed the normal linear variation in 1/H above about 8 kilo-oersteds and 
extrapolated to a value at (1/H=0) 18 per cent. higher than at H=10-6 kilo- 
oersteds. Comparison of the torque curves of specimen 3 at 6-5 and 10-6 kilo- 
oersteds showed that it is only the sin 20 term which increases with increasing 
field. The sin 46 and sin 66 terms decrease in magnitude. The decrease in 7’, 
is so strong that this term extrapolates to zero at (1/H=0) within the limits of 
experimental error. The extrapolated ratio 7,/T, is approximately —0-07. 
However, the striking, qualitative difference between the torque curves of the 
iron and nickel specimens is not influenced by the magnitude of the field in which 
the measurements are made. 


There is another experimentally variable parameter in the particular case 
of nickel, in which the first magnetocrystalline anisotropy constant is negative 
at room temperature but becomes positive (iron-like) above about 100 °C. 
Specimens 7 and 8 were made by mixing weighed samples of nickel powder with 
pyrex powder and either benzoic acid or ‘‘ Armo wax’ (Armour Chemical 
Division) to form a thick paste at about 150 °C; the specimens were then allowed 
to cool in beakers in the 10-6 kilo-oersted field, so that they solidified well above 
100 °C (see Table 1). The anisotropies of the specimens of nickel powder set in 
this way were more than 1000 times stronger than for the same grade of powder 
set in plaster at room temperature. 


Discussion 


The shape of the torque curves can be explained almost entirely in terms of a 
stringing together of grains. It is surprising therefore that of the specimens 
which had been set at room temperature only those containing iron showed this 
effect. A correlation with magnetocrystalline anisotropy is indicated by the 
torque curves of specimens containing nickel powder which had been set above 
100 °C (7 and 8) and which showed very strong anisotropy similar to that of 
the iron specimens. Above about 100 °C the first magnetocrystalline anisotropy 
constant becomes positive (iron-like) (Honda, Masumoto, and Shirakawa 1935 ; 
Bozorth 1951). 

A physical difference between the iron and nickel powders must be noted. 
Carbonyl iron powders are polycrystalline, but when nickel powder is produced 
‘the particles tend to be single crystals and to possess an idiomorphic shape ”’ 
(Pfeil 1947). Microscopic examination of the magnetite grains used in specimen 9 
showed that these also had irregular but characteristic Shapes indicative of 
single crystals. The conclusion may therefore be advanced that polycrystalline 
magnetic particles (iron, specimens 2, 3, 4) or single crystals with positive magneto- 


SHORT COMMUNICATIONS 201 


crystalline anisotropy (nickel set above 100 °C, specimens 7 , 8) tend to string 
together along lines of force when exposed to a magnetic field. Single crystals 
with negative magnetocrystalline anisotropy (nickel below 100 °C, specimens 
5, 6, and magnetite, specimen 9) do not show this effect. No mechanical 
explanation for this difference in behaviour can be offered at present. 


The torque curves can also be examined for details indicative of micro- 
crystalline alignment. For material with positive magnetocrystalline anisotropy 
this would be most apparent as an increase in the sin 40 term in the torque curve 
by the addition of a term having the form of equation (2) to the curve due to the 
stringing effect. This means a decrease in magnitude of the negative coefficient 
T, relative to the magnitude of 7,. Experimental uncertainty in determining 
the ratio 7,,/T, (see Table 1) does not permit any distinction to be made between 
specimens 2, 3, and 4, so that carbonyl iron powders of grades MCP and MF are 
evidently both polycrystalline, as mentioned above. 


The occurrence of a magnetocrystalline term in the torque curves of 
specimens 7 and 8 appears possible, particularly in the latter. At high temper- 
atures the nickel powders would experience [100] crystal alignment in the field, 
so that with this alignment frozen in at room temperature the addition to the 
torque curve of a term having the form of equation (2) with negative K, will 
increase the magnitude of the negative 7, term relative to T7,. The increase in 
(—T,/T,) which was observed for specimen 8 is well outside experimental error 
and is consistent with crystal alignment in the nickel powders at high 
temperatures. 

The difference between specimens 5 and 6 is evidently the result of a real 
physical difference between the nickel powders used. The Sherritt-Gordon 
grade C50 powder gave no measurable anisotropy at all, whereas exactly the 
same treatment of a sample of Mond carbonyl grade A powder resulted in a 
specimen of very small but readily measurable anisotropy. If the foregoing 
conclusion that anisotropy will only appear if the grains are polycrystalline 
or have positive magnetocrystalline anisotropy is correct, then it appears that the 
Sherritt-Gordon powder is entirely composed of single crystal grains, but that the 
Mond powder contains a small but noticeable fraction of polycrystalline grains. 


Powders for these experiments were provided by the Mond Nickel Company, 
London, and Sherritt-Gordon Mines Ltd., Alberta, Canada. 
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ELECTRON BACKGROUND IN NUCLEAR EMULSIONS* 
By V. D. Hopper} and JEAN E. LAByt 


Since 1951 supplies of Ilford G5 electron sensitive nuclear emulsions have 
been flown by air to Australia. In 1954 it was observed that some batches, 
when developed soon after arrival, showed an excess electron background which 
made them unsuitable for the cosmic ray work for which they were imported. 
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Fig. 1.—Dosage in milliréntgens received by electron sensitive emulsions arriving by air 
frona England (A), and arriving by ship (S) over the period 1955 to 1959 (June). 


It was then decided that a test development was required for all emulsions 
before they could be Safely used in the high altitude studies. It was also con- 
Sidered necessary to import samples by ship to determine whether emulsions 
would be more suitable if imported by this means. 


* Manuscript received January 5, 1960. 
{ Physics Department, University of Melbourne. 
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Results of the dosage received by samples arriving throughout the period 
January 1955-June 1959 are given in Figure 1, where plates sent by air are marked 
A and those by sea 8S. The electron background dosage in milliréntgens 
calibrated as described earlier} is plotted against the approximate date of manu- 
facture. Correlation of the heavy dosage Samples with radioactive rainfall in 
France was noted earlier.* 

The emulsions consisted of 600 micron thick Ilford G5 emulsion on glass 
backing. The size of the plates was 3 by 2 in., and these were packed in cardboard 
packets holding 6 or 12 plates. The packets were in turn protected with a layer 
of wood shavings placed inside wooden boxes. These boxes were labelled photo- 
graphic material and the label specified that the packet should be kept 10 ft 
from any radioactive material. Comments regarding radiation due to instrument 
panels etc. are made in the previous communication on this subject. The same 
method of packing has been made throughout the duration of this test. 
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Fig. 2.—Histogram showing the frequency of dosage: (above) plates 
arriving by ship, (below) plates arriving by air. 


It has been the standard practice of Ilford Ltd. to dispatch plates within 
about 3 days of manufacture and to mark the manufacture date on the box 
carrying the plates. | 

Figure 2 shows the histograms for plates arriving by sea and by air and 
clearly shows that air-flown samples are more prone to excess radiation damage 
than those arriving by sea. (The mean air transport time was 38hr at 
~ 10,000-20,000 ft, climbing and descending time about 8 hr.) 

Only two samples out of 60 of those arriving by sea showed an exposure 
exceeding 25 mr. One sample was manufactured in September 1955 and the 
corresponding sample which arrived by air also showed excessive exposure. 


* Hopper, V. D., and Lasy, Juan EB. (1955).—Aust. J. Phys. 8: 557. 
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The other sample, manufactured on March 2, 1959, showed a very high exposure 
(170 mr), whereas the corresponding sample arriving by air showed low exposure 
(13 mr). The reason for this difference has not been determined. 

Results show that it is advisable to import nuclear emulsion plates by sea 
rather than by air, although there is still some chance of excessive exposure. 
The use of the electron sensitive emulsions in radiation measurements is also 
shown to be useful for determination of low dosages, particularly as the electron 
tracks can be distinguished from darkening produced by other means such as 
exposure to light and chemical fogging due to excess temperatures. Further 
tests are being carried out with the method on aircraft travelling on different 
routes in the southern and northern hemispheres. 


RELATIVE INTENSITY OF THE 17-2 AND 14-3 MEV GAMMA 
RAYS FROM THE “Li(p,y)*Be REACTION* 


By B. MAINSBRIDGET 


Radiation from proton capture in 7Li is known to consist of two principal 
components of energy (17:2+ jH,) and (14-3+{H,) MeV, corresponding to 
transitions to the ground state and first excited state of §Be respectively (Walker 
and McDaniel 1948). Resonances in the reaction are known to exist at 
E,=441 keV, 1-03 and 2-1 MeV (Bonner and Evans 1948; Kraus 1954; Price 
1954) and the relative intensity of the two y-rays is known to vary in the 
neighbourhood of the 441 keV resonance (Campbell 1956). It is not known if 
the intensity ratio varies in the region of the 1030 keV resonance and this experi- 
ment was designed to repeat the measurements of Campbell and extend the 
investigation to the higher resonance. 


Haxpervmental 


The detector was a 5in. diameter by 4 in. long Nal(Tl) crystal mounted, 
with the cylindrical axis vertical, at 90° to the proton beam and 9-5 in. from 
a 7s in. diameter natural lithium target. The y-radiation was collimated into 
the crystal centre by a conical hole through a lead block 3 in. thick. The col- 
limator was tapered so that the detector observed, through a solid angle of 
0-074 steradians, a + in. diameter spot at the target. A 5 in. diameter Dumont 
6364 photomultiplier was mounted on the crystal and the output was amplified 
and fed to a 100 channel pulse-height analyser, biased at 6 MeV y-ray energy. 


* Manuscript received March 7, 1960. 
+ Research School of Physical Sciences, Australian National University, Canberra. 
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The targets, prepared by the evaporation of natural lithium metal onto 
silver-plated brass backings, were mounted on a water-cooled holder and inclined 
at 20° to the proton beam. The deposits were as thin as experimental conditions 
would ailow and were less than 5 keV thick for measurements between E,=430 
and 450 keV. Up to #,=600 keV, 5-10 keV thick targets were used and in the 
remaining energy regions they were 20 keV thick except for the measurements 
below 300 keV, where the low yield required targets 60 keV thick. 
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Fig. 1.—Pulse-height distributions of the 7Li(p,y)®Be radiation at By, =440, 620, and 
1030 keV. 


Beam currents of 10-15 pA were used and the charge was measured on a 


current integrator. In the region of the 441 keV resonance, the current was 


reduced to limit counting rates to 100 p.p.s. 
The bias and gain of the pulse-height analyser were fixed during the measure- 
ment so that the peak of the 17-6 MeV y-ray always occurred in the same chee 
of the display. The cosmic ray background was measured by running the 
electronics for 2-3 hr with the accelerator switched off. 
Figure 1 shows typical spectra obtained in the measurement after ates 
for cosmic ray background and normalizing to a common peak height at 17-6 MeV. 
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Pulse-height Distribution of the 17-6 MeV y-Ray : 

The pulse-height distributions of the two y-rays cannot be obtained directly 
from the spectra of Figure 1 because of the overlap in the two components. 
The separation is complicated by the 2 MeV width of the 14-8 MeV y-ray. 

The 17-6MeV component is monochromatic and investigations of the 
response of 5in. diameter by 4 in. long Nal(Tl) crystals to monochromatic 
radiation between 10 and 20 MeV suggest there is little change in the pulse-height 
distributions in this energy range (Kockum and Starfelt 1959). A comparison 
of the pulse-height distributions of the 16-6 and 20-3 MeV y-rays from proton 
capture in “B and °H respectively, measured with the apparatus of the present 
investigation, showed the line shapes to be identical and if was assumed that the 
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Fig. 2.—-Variation of the intensity ratio with proton energy, measured at 
90° to the proton beam. 


pulse-height distribution of the 20-3 MeV y-ray could be used to represent 
that of the 17-6 MeV y-ray. This was obtained by replacing the lithium target 
with one of tritium absorbed in zirconium and bombarding with 850 keV protons. 
The gain of the photomultiplier was altered so that the peak of the 20-3 MeV 


y-ray coincided with that of the 17-6 MeV component and the line shape is 
included in Figure 1. 


Line Shape below Hy,=8-5 MeV 

The pulse-height distributions do not indicate the line shape below 
Hy,=8:5 MeV. In this region, contributions from the 9H (p,a,y)!®O reactions 
are mixed with low energy y-rays from scattered radiation. Many of these 
events arise from scattering material outside the detector and some are due to 
the escape of a fraction of the incident y-ray energy from the detector. The 
actual line shape between 0 and 8-5 MeV depends on the escape loss and is 
partly determined by the dimensions of the collimator. It was inferred from 
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the investigations of Kockum and Starfelt (1959) that the line shape for the 
geometry of the present experiment was between the limits in the extrapolations 
to zero energy in Figure 1. 


Calculation of the Intensity Ratio 

Each spectrum was summed between #,=8-5 and 19-5 MeV and the cosmic 
ray background was subtracted. The assumed pulse-height distribution of the 
17-6 MeV y-ray was normalized to fit each spectrum and the intensity of the 
14-8 MeV component was obtained by subtraction. Small corrections were 
made for the calculated efficiency of the detector to 14-8 and 17-6 MeV y-rays, 
82-9 and 84-5 per cent. respectively. The ratio of the intensity of the two 
components I,,.4/Z,4., is shown as a function of proton energy in Figure 2. The 
excitation function was derived from the ratio of the total number of events 
between 8:5 and 19-5 MeV and the current integrator counts. Four targets 
of varying thickness were used in the measurements and, since the yield from 
each target at the 441 keV resonance has been normalized to a common value, 
the excitation function is only approximate. 


Absolute Value of the Ratio 

The experimental points in Figure 2 are based on the pulse-height 
distributions above E,—8-5 MeV. The errors in the absolute values of the 
ratio, due to the uncertainty in the line shapes below this energy, have been 
calculated by considering the extrapolations of Figure 1. The absolute value 
lies between the experimental curve and the broken curve of Figure 2. 


TABLE | 

MEASUREMENTS OF THE INTENSITY RATIO AT THE 441 keV RESONANCE 

Ky, aot Angle Ly 7-6/Lia-8 Reference 
Thickness 
450 10-20 keV (Ue, Chie 2+0-2 | Devons and Lindsey (1950) 
460 150 keV 0° 2 Walker and McDaniel (1948) 
500 thick Os) Sie We 1-7+0:2 Stearns and McDaniel (1951) 
44] 5 90° 2-30+0-04 | Present work 
Discussion 


The points at H,—480 and 500 keV show higher values of the ratio than 
neighbouring points and this effect is attributed to non-uniformity of the target 
deposit where a small portion may be thick enough to contribute some resonance 
radiation with a larger cross section and a higher value of the ratio. In this 
case, the target was discarded and the measurement was repeated. The results 
are in general agreement with the measurements by Campbell (1956). The 
improved resolution of the larger detector used in the present investigation limits 
the errors due to uncertainty in the line shape of the two y-rays and allows the 
ratio to be measured directly without reference to measurements by other workers. 
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At the 441 keV resonance, the value of 2-30-04 is higher than values 
reported by other workers using pair spectrometers (Table 1), but, since those 
measurements were made with thicker targets, a small contribution from the 
non-resonant radiation would tend to reduce the value. 

Above 800 keV, the ratio remained constant at 0-54-+0-08 and no departures 
from this value were observed at the 1030 keV resonance. 

Measurements of the angular distributions of the two y-rays have shown the 
intensity ratio to be in addition a function of the angle of observation above 
and below the 441 keV resonance and in the neighbourhood of the 1030 keV 
resonance and is the subject of further investigation. 


The author wishes to thank Professor E. W. Titterton for his interest in the 
programme and the Australian Atomic Energy Commission for a grant which 
made the work possible. 
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